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Abstract 

We discuss in this set of lectures the structure of proton/neutron as re- 
vealed through a study of form- factors. This is followed by a discussion of the 
structure functions of proton/nuclear targets as measured in the deep inelas- 
tic scattering (DIS) of leptons off these targets. We discuss the parton model 
in DIS as well as outline the usage of parton model in processes other than 
the DIS. We then go on to dicsuss the EMC effect: the nuclear dependence 
of the structure functions. After a discussion of different models of the EMC 
effect we end by pointing out the possibility of distinguishing between these 
different models by studying the correlation between the j4-dependence of 
different hard processes and the EMC effect. 



*Lectures given at the S.E.R.C. school in Nuclear Physics, Goa University, Feb. 1993 



Introduction 

Structure of hadrons, the strongly interacting particles, has been the subject of the- 
oretical and experimental investigations over the period of past four decades and has 
played a crucial role in the birth of the subject of particle physics. Traditionally, in- 
formation about the constituents of a system has come from two sources: one is the 
study of static properties such as mass, magnetic moments, spin, parity, etc., and 
the other is scattering experiments. The famous Rutherford scattering experiment 
is the prototype of the latter. In the case of hadrons, the clues to quark structure 
came in the form of the quark model put forward by Gell-Mann |I| , but the final con- 
firmation came from the Deep Inelastic Scattering (DIS) experiments where the 
observed phenomenon of 'scaling' was explained |^ in terms of 'partons' (pointlike 
constituents of proton) [Q. The discovery of 'asymptotic freedom' of 'Quantum 
Chromodynamics' (QCD) ^j: the gauge theory of quarks and gluons, coupled with 
the operator product expansion gave an explanation, in the context of perturbative 
QCD(pQCD), why the parton model works so well[0]. By now QCD is the accepted 
theory of strong interactions. In the set of these six lectures an effort will be made 
to review the phenomenological understanding of DIS and how its study led to the 
parton model, along with some aspects of parton structure of nuclei. We begin with 
an introduction as to how one probes the structure of a nucleus via electromagnetic 
interactions and explain the concept of the electric and magnetic form factors of 
the proton/neutron. Following this, we will discuss the DIS kinematics and intro- 
duce the idea of two structure functions in terms of which the electromagnetic DIS 
cross-section can be parametrised. This will be followed by a discussion of scaling 
of the structure functions and parton model. The scaling violations predicted by 
QCD will be alluded to very briefly. In the case of the DIS processes mediated 
by the electroweak gauge bosons, W^/Z'^, three independent structure functions 
are required to parametrise the DIS cross-section. We introduce these and briefly 
discuss the important role played by neutrino induced DIS processes in providing 
conclusive evidence for the parton model. The momentum distribution of partons 
measured in both the electromagnetic and the weak DIS processes is found to be 
affected by the nuclear environment. This is the so called EMC effect 0, |^. This 



will be discussed in the last lecture. Various models [y, |T0| have been put forward 
as a theoretical explanation of this effect. Some of the models will be summarized 
along with a discussion of possible experimental tests to distinguish among these 
models [p^ . 

1 Electromagnetic Structure of Hadrons 

Historically, hadrons were classifled as particles with large masses (m ~ 1 GeV) 
which participate in strong interactions. The wide variety of these hadrons {p, n, A, 
S'^, S^, S^ , TT^, 7r°, K^, K^, K^, ■ ■ ■) already indicated that these are not elementary. 
Further, for an elementary, spin ^ particle, Dirac theory predicts gyromagnetic ratio 
(yf to be 2. Even the deviation of g from 2 can be computed in perturbation theory. 



For example, for an e , theoretical computations |]TT| predict 

^ ^ -0.328478966 - + 1.1765(13) - - .8(2.5) - 



2 ) e 27r \7r/ \tx ) \tx ) 

= 0.001159652460(192) (1.1) 

where a is the fine structure constant given by e^/An.f^ The experimentally measured 
is 

'^-2yP*- 



value 12 



-) =0.001159652193(10) (1.2) 

/ e 



2 
In eqs. ( |1 . 1| ) and ( |1.2| ) the numbers in the bracket indicate the theoretical and 



experimental uncertainty respectively. The excellent agreement between the theo- 
retical prediction and the experimentally measured value indicates that the electron 



is indeed elementary. However, for a proton and neutron, the measured values ||12 
of magnetic moments are 

fip = 2.79 fiN = 2.79 



2Mp 



/i„ = -1.91 /i^ = -1.91^ ■ (1.3) 

These correspond to values of 1.79 fiN and —1.91 fiN for anomalous magnetic mo- 
ments of the proton and neutron respectively. This clearly indicates that they are 
not elementary particles i.e. they have spatial extension. They have a structure. 
But this is, strictly speaking, an energy dependent statement. At low energies, the 
scattering experiments are not able to 'resolve' this structure. After all, Rutherford 



concluded from his scattering experiments |jT3l with a particles that all the positive 
charge in an atom is concentrated in a point nucleus and the electrons occupy the 
rest of the atomic volume. But what this conclusion really meant was that the 

o 

nuclear radius, r^, is much smaller than atomic sizes (~ O (few A f» 10~^° m)) 



One had to perform the scattering experiments |T^ with higher energy electrons 



(~ 100 — 500 MeV) to reveal the structure of nuclei at a distance scale ~ few fm 
= 10-^^ m. 

1.1 Effect of structure of the scattering centers on scattering 
amplitudes 

To understand this let us consider a general scattering process 

e~ + A^e'+A. (1.4) 

The process is represented digramatically in fig. |1]. Let Pi and P2Q denote the four 
momentum of the target nucleus and the incoming electron respectively and P3, P4 

^I will use units h — c — I throughout these lectures. 

■''I use Pauli metric such that P^ = ~"if . Details of notation are given in Appendix A. 



be the four momenta of the e and nucleus in the final state, respectively. Let 

Pi = {0,zMa), P2 = (Po,^^o), 

Ps = {p,iE), P,= {p',iW) . (1.5) 

As shown in the figure, the electromagnetic scattering takes place via the exchange 
of a photon. Application of four-momentum conservation at each vertex implies 
that the four momentum of the photon is given by 

q = P^-P^ = P^-P, . (1.6) 



In the notation of eq. (|1.5| ) we have, 

q = (g, iqo) = {f - P, i{Eo - E)) 

= {p\ i{W-Ma)). (1.7) 

Using eqs. (|1.5|) and ( p..6|) , we get for the invariant mass of the photon, 

g2 = (P2 - P3)2 = -2ml - 2IP2I \pz\ cose + 2E2E3 . (1.8) 

Neglecting the electron mass rrie (which also implies \p2\ = Eq = \po\ = po and 
IpsI = E = p) we get, 

q = 2ppq{1 — cos 9) = AppQ sin - ' (1.9) 

Thus for a scattering process one always has g^ > 0. This means that the exchanged 
photon is not a real photon but a virtual one. 

It can be shown, on very general grounds, that the net effect of the presence of 
a structure in the scattering center on the scattering amplitude, calculated from a 
diagram such as shown in fig. |1], is to multiply it by a form factor F{q^). This form 
factor is given by, 

F{q')= le^'-y f{y)<fy (1.10) 



where f{y) is the distribution function describing the target. In the limit of recoil- 
less, elastic scattering {W ~ Ma) eq. ( |1.6| ) gives q = {q,iO). In this case the 
four-fourier transform reduces to the spatial Fourier transform. 

1.2 Relation of spatial charge distribution to form factor 

Let us consider scattering of a spinless electron from a spinless charge distribution 
and derive an expression for the form factor in this case. Consider an e~ of energy 
£"0 and momentum po incident on a target nucleus of charge Z\e\ at rest. The four 
momenta of various particles are then given by eq. ( p..5|) . For a point nucleus, the 
classical Rutherford scattering formula for the differential cross-section^ is given by 



da 

In 



Aptsm'l 4p^sin^f ^ ^ 

Rutherford " ^ " ^ 



^Note that this formula can be derived classicaUy as weU as in non-relativistic quantum me- 
chanics using Born approximation for the scattering of an e~ in the screened, Coulomb field of a 
nucleus. 



where Q is the sohd angle and dQ = smOdOdcj), 9 being the scattering angle. For 
a spinless, relativistic e~ scattered by the screened, Coulomb potential of a point 
nucleus, this expression is, 

"" ''"' • (1.12) 



dn AEl sin^ f 

If we assume the initial electron to be incident along the z axis, the kinematics is 
as depicted in fig. ^ Now using eqs. ( |1.5| ) and ( |1.6| ), we see that 

Po=P + p' , (1.13a) 

Ma + Eo = E + W . (1.13b) 

If Ei and Ef are used to indicate total initial and final state energies, we have 

Ei = Ma + Eo = Ma + Po = Ef = W + E = W + p . (1.14) 

Let us suppose that the nuclear charge Z\e\ is distributed with a distribution function 
p{R), normalised to unity by, 

'^ p{R)d^R= 1 . (1.15) 

Hence the total electric charge in a volume element d^R is given by Z\e\p{R)d^R. 
The scattering cross-section is 

W 
da = — , (1.16) 

V 

where W is the transition probability and v is the velocity of the incident e~ w.r.t. 
the center (which for relativistic electrons is 1 in our units) .|^ The transition prob- 
ability W is calculated using Fermi's golden rule 

W = 2npf\Mif\'^ , (1.17) 

where pf is the density of states in the final state and Aiif is the scattering ampli- 
tude, computed in Born approximation for the screened. Coulomb nuclear potential. 
For the two particle final states, the density of final states pf is given by (with h = 1), 

dN 2 JO dp ., ,„. 

and Aiif is given in the Born approximation as 

Mif = f ip}ir)V{r)ipi{f)d^r . (1.19) 

dp/dEf in eq. ( |1.18| ) can be computed using eq. ( |1.13a| ) as follows. Eq. ( |1.13aj ) 
gives 



Ef = p + W = ylp'l'^ + rriA +p 



Ma + Pq + p"^ — 2ppo cos 9 + p 



"Note here that I am normahsing with one particle per unit volume. 



This in turn gives, 



W 



dp 

dEf Ef — po cos 9 



;i.2o) 



The energy and three momentum conservation of eqs. (|1.13a|) , (|1.13b| ) gives 



Ef = Po + Ma = p + W = p + y p2 + p2 _ 2pp^ cos 9 + m\ . 
This gives us, 



— Ma = Ef — pq cos 9 . 
p 



This can also be rewritten as 



P 



Po 1 + 2po/Ma sin 



2 e 



Combining eqs. (|1.2(]| ) and (|1.22|) we have 

dp W p 



dEf Po Ma 



(1.21) 



:i.22) 



:i.23) 



;i.24) 



To calculate Aitf in eq. ( |1.19| ) we need V{r). If we take into account the 
screening of the nuclear charge due to atomic electrons, the potential V{r) felt by 
the electron at f, due to the nuclear charge distribution shown in fig. S will be. 



V{r) 



Ze^ r p{R)d^R 
4n J \f-R\ ^ 



-|r-_R|/a 



1.25) 



where a is a damping factor ~ O (atomic radius) which arises from the screening 
of the nuclear charge by the atomic electrons and is called the screening radius. 
Clearly a ^ R. If one uses the plane wave approximation for the wave functions of 
the incident and scattered electrons, with momenta po and p respectively, we have 



rAr) 



'tp-r 



Ipiif) 



ipo-r 



{h = c = i) . 



1.26) 



This then 


gives for M 


f ' 




M.f = 


-Ze^ 






-Ze^ 



\f-R\ 



e'^-^rfV 



d^R 



iq{r-Ry\r^R\/a 

e''^-^p{R)d'^R I ^ d^R 

If— R\ 



4:71 

Choosing the z axis along g and s = f—R, we get 



M 



if 



Atx 



e'^-^p{R)d^R / 2 7C sds dcosa e*'?'™'^" e""/" 



:i.27) 



The first integral is a function only of Q^ = \q\^. Then we have 



^^f = -Fm^ r.. , .,.. ' (1-28) 

where 



47r g2 + 1/^2 ' 



F(Q2) = / e'^-^p{R)SR . (1.29) 



F{Q'^) is called the form factor. Note here that the screening factor e~^/" in the 
integrand is essential to make this integral convergent. Q^ = |^^ = PQ+p'^—2ppQ cos 6 
is ~ (9(MeV^) for the energies under consideration whereas a ~ O{10~^cm) ; z.e. 
l/a^ ~ 4(KeV^). | Hence for Eq^O (MeV) , l/a^ < g^ Therefore the expression 
for M.if can be approximated as 






M.f^—^^-^. (1.30) 



Substituting eqs. ([Ool) , ([T^OD and ( ITTsD in ([lTtD , we get 



,,.2.P!f!^_Jl_.f££!yM!|f,g'.)|. . (1.31, 

(27r)3Ej-poCOS^ V^TT y Q4 I v^ ^1 y ) 

Using eq. (|1.23|) we get finally, 



,2\ 2 I T7/,o2\|2 



da W p Zey \FiQ^ 

Since the nuclear mass Ma ~ O (GeV) ~ (9(1000 MeV), for electron beam energies 
of ~ C (MeV) , W ~ Ma ■ Hence eq. (|1.13bD tells us p ^ po and hence Q"^ ~ 
4pg sin^ 6/2 = AEq sin^ 6/2. Using all this the expression for the differential cross- 
section becomes 

wm? ■ (1-33) 



da Z'^c? 



dn 4E^ sin^ ^ 



Comparing eqs. (|1.33|) and (|1.12|) we see that the net effect of the distribution of 



the nuclear charge over a spatial volume according to a distribution function p{R) 
is to multiply the cross-section by the square of the form factor |F(Q^)p. F{Q^) 
is nothing but the three dimensional Fourier transform of the charge distribution 
p{R). If |F(Q2)| = 1 for all Q^ then eq. ( [03D reduces to eq. (|T:T2|) which is 



the Rutherford scattering cross-section formula for the scattering of a relativistic 
electron from a spinless point nucleus; i.e., 

\ / charge distn. \ / point nucleus 

II Use he — 197 MeV fm to arrive at this result. 



Table 1: The form factors and the Root Mean Square (r.m.s.) radii (i?^)^/^ for 
different charge distributions. 



p{R) 


Fm 


(i?2)V2 




6^{R) 


1 


0.0 






1 


V6 
m 


Monopole 


(1+4) 


Stt 


1 


VV2 
m 


Dipole 


(1+4)2 



where the form factor is given by eq. ( |1.29D . For a spherically symmetric charge 
distribution, choosing the z axis along q, the equation for F{Q'^) becomes 



F{Q' 



'^L ^(^) QR 



'"^^^RHR 



;i.35) 



The charge distribution p{R) is given in terms of F{Q^) by the inverse Fourier 
transform, 

P(R) = ^2 J FiQ'f^Q'dQ . (1.36) 

At values of Q = l^l such that QR <C 1 over all the region where p{R) is appreciable, 
one can expand sin QR/QR in the integrand in eq. ( p..35| ) and we get 



fOC I )^ POO 

F{Q^) = Atc / R^p{R)dR - ^ / A'KR^p{R)R^dR + OiQ^R") . 
Jo 6 JO 

Hence for a charge distribution normalised as given by eq. ( |1.15|) we have, 

Q'{R') 



1.37) 



F{Q') = 1 



6 



:i.38) 



Thus a measurement of F{Q'^) at small Q"^ such that Q'^{R^) <^ 1 gives us {R^) for 
different charge distributions. Table |l| summarises the form factors F{Q'^) and (i?^) 
for different charge distributions where these can be calculated analytically. As can 
be seen, all of them satisfy eq. ( |1.38| ). 

Recall now from eq. (|1.7| ) that the square of the four-momentum transfer q^ is 
given by 



\P 



•/|2 



{W - Ma] 



-2MX + 2WMa . 



;i.39) 



This gives us 



q^ 



2Ml 



+ 1 



W 



;i.4o) 



Further for W ~ Ma [i-e. q^ -C 2m^; small nuclear recoil), we also have q^ ~ Q^. 
This means that in this approximation we can replace Q^ by q^ in all the earlier 
equations. 

If we repeat the analysis above for a point e~ incident on a point spinless charge 
Z\e\ but keeping the effect of e~ spin we get (see problem 4, Appendix B) 




Z^a^cos^l 



point e 
point, spinless nucleus 



4E2sin^f 



1 + fe-n^f 



'1.411 



Neglecting the recoil, for relativistic e -nucleus scattering this becomes, 

_ Z^a" cos2 f _ AZ^a^ cos^ f E^ 




Mott. 



4^.2 sin^ 



:i.42) 



The factor of cos^ | in the numerator in eq. (|1.42| ) indicates the impossibility of 180° 
scattering for longitudinally polarised, spin | electrons (see problem 3, Appendix 
B for a further discussion of this point). Once again, repeating the exercise for a 
nuclear charge distribution, we get 



dn 



spin Tye 
spinless, charge distn. 




\F{q 



2m2 



;i.43) 



Mott 



where F{q'^) is the form factor given by eq. ( |1.29| ). Hence for a spinless nucleus the 
form factor can be measured by the ratio: 



2\|2 



\n<i')\ 



(da/dQ) 



eA^eA 



{da/dn)uott 



:i.44) 



The charge distribution of the nucleus is then given by the inverse Fourier transform 
of F(g^). Hofstadter studied fl^ the e~-nucleus scattering for ii^e-~600 MeV. Hence 
the electrons were relativistic. Since M4 ~ 1000 A MeV, the slow recoil assumption 
is also justified in this case. Hence, the nuclear form factors could be determined 
in this case using eq. ( |1.44| ). Thus a deviation of the angular dependence of the 
cross-section from the one expected for a point target, measures the form-factor. 
Therefore to get information about a possible structure of the scattering center, we 
must know the theoretical predictions for the results expected for a point target. 
These, along with the experimental information on the form-factor, for the case of 
proton/neutron will be discussed in the next lecture. It must be noted here that 
the interpretation of the form-factor as a three dimensional Fourier transform of the 
charge distribution is strictly true only in the case of recoilless scattering. At higher 
energies, when this assumption is not justified, the form-factor F{q^) can be looked 
upon as the ratio of the scattering amplitudes from an extended and point target. 



2 Form Factors of Proton and Neutron 

2.1 e j9 -^ e^p for pointlike "Dirac" proton 



The discussion in the last section and eq. ( |1.3|) in particular indicates that the 
proton and neutron are not pointlike and elementary. It also makes it clear that 
one can obtain information about the charge distribution in a proton (neutron) by 
studying the electromagnetic scattering process 

e~ + p (n) —>■ e^ + p (n) . (2.1) 



Eq. (|1.44|) indicates that to be able to do this it is essential to know the expected 
cross-section da/dfl for a point proton. In the discussion below we will try to 
indicate how such a calculation is done. 

At the energies which we are considering the proton and neutron also need to be 
treated relativistically. For a Dirac (i.e. point) e~{p), the electromagnetic current 
is given by 

= ,^^lM^(9l^^^m-,^m^^(^a..^^) , (2.2) 



2meiMp)\dx/ ^dxj 2meiMp)dx,^^ ""^^ 



where 7^ are the usual Dirac matrices and a^^i, = ^[7/i,7i^]; Qeip) is the charge of the 
e~{p) and me{Mp) is the mass of the e~{p). The fermion index (e~ or p) on the spinor 
is suppressed. It can be shown easily that this expression for the electromagnetic 
current means (^e = 2, /ig = 1 B.m., Qp = 1 and /ip = 1 n.m. To calculate da /(El one 
needs the matrix element for the scattering process e~p —>■ e~p. This is given by 

M ~ ^J^J^ , (2.3) 

where g^ is the square of the four momentum transfer in this process. The Feyman 
diagram for the scattering process and the four momentum assignments of various 
particles involved in it are shown in fig. |^. The kinematics of the process is the 
same as given by eqs. ( p..5|) - (|1.9| ), (|1.134 |l-13b|) and (|1.22|) - ( p..24|) , after replacing 



Ma by Mp. Again the square of invariant mass of exchanged photon g 7^ 0, and 
hence it is virtual. 

To calculate da/dQ now one has to use the rules of quantum field theory. The 
cross-section for a process A + B—>-C + D is given by 

1 1 (27,)'S*(Pa + Pb- Pc - Pd) iPpc <PPD 



where Pa, Pb, Pc and Pd are the four momenta of the four particles involved and 
|A^p is the square of the matrix element averaged over initial state spins and 



summed over the final state spins. The matrix element Ai itself for the scatter- 
ing process is given by eq. ([^.31). Using J^^^^ given by eq. ( ^.21 ) and the trace rules 
given in Appendix A , we get for a point proton (see problem 4, Appendix B) , 



da 



>) 



a^ cos^ I 



Point e-/p 4^2 Sin^ f 1 + ^ siu^ f 



1 + 



2M^ 



2^ 

tan - 




p_ 

Mott Pq 
(Z=l) 



1 + 



2Mp2 



tan 



(2.5) 



The first term is just the Mott electrostatic scattering cross-section of eq. ( |1.42|) 
with the nuclear charge Z = 1, the factor of p/po is the recoil factor which goes to 
unity for recoil- less scattering (recall eq. ( |1.13b|) ) and the last factor is due to the 
magnetic moment of the proton and electron which reduces to unity if proton were 
to be replaced by a spinless point charge. Thus for a scalar, point proton we can 
write (see problem 4, Appendix B), 



da 

dn' 



,e p 



>) 





a^ cos^ 1 1 


[ 2MI 


int e~/p 


4^2 sin" f 1 + g sin2 | 




(da\ p 




\dQ, J Mott Po 
(Z=l) 





tan 



e 



(2.6) 



Note also that we are no longer restricted to ly ~ Mp and hence q^ ^\q\^ = Q"^ 



2.2 Effect of anomalous magnetic moment of proton 

The above discussion assumes that the proton is a Dirac particle with gp = 2. In 
reality, even if we were to neglect the possibility that the proton is an extended 
object, the anomalous magnetic moment of the proton cannot be neglected, i.e., eq. 
( p.2|) cannot give the electromagnetic current for the physical proton even if were to 
be point like. The anomalous magnetic moment gives an additional contribution to 
the electromagnetic current of a proton and it can be written asQ 



' TP " 



_iqp_ 
2M, 



P L 



dXn ^ 






/^J 



qp 

2M„ 



[1 + Kp) 



^ T / 



dx 



(2.7) 



Note here that the effects of a possible spatial structure of the proton are not taken 
into account apart from the nonzero value of the anomalous magnetic moment. 
The interaction Hamiltonian ?i*"* is given by ~ "J^ " A^ and hence the scattering 
amplitude will be determined by this. Recall also ipix) ~ u{P) exp{iP ■ x). Since Pi 
and P4 are the four momenta of the initial and final state respectively, the expression 
for the current of eq. (|2.7| ), in momentum space, becomes 



" TP 



J. 



Dirac 



+ 



IQp _ 



2Mr 



Up{Pi)a^^Up{Pi)q^ 



(2i 



**see for example Advanced Quantum Mechanics, J.J. Sakurai, (Addison Wesley, Reading, Mass. 
USA) 



where J^*'''^'= represents the current of eq. ( p.2| ), in momentum space, given by 
Jf ™^ = iqr>u,{PA)l,u,{P,) 

- "^^ {Pa^ + Pi^.)up{Pi)up{P{) + ^q,Up{Pi)a,^Up{Pi). (2.9) 



2M/ "'^ -.^,^^v-.,^^v-., 2Mp 

This is called the Gordon decomposition of the electromagnetic current of a spin 
1/2 particle. 

2.3 Effect of the spatial extension of proton 

The discussion upto now has neglected a possible spatial extension for the proton. 
However, the large value of Kp given by eq. ( |1.3| ) implies that the proton is not a 
spin 1/2, point particle. Our discussions on nuclear form factor indicate that the 
proton structure will be reflected in a multiplicative form factor for the scattering 
amplitude. Following this, we can try to write down the most general expression 
for the electromagnetic current of a spin 1/2 object by constructing a four vector 
using the various bilinear covariants made up of '^/'^(-Ri) and ippiPi) along with the 
independent four momenta available in the scattering process of eq. ( p.l|) . The 
requirement that parity be conserved, i.e., the electromagnetic current be a vector 
and not a pseudovector, rules out all the bilinears containing 75. The possible 4- 
vectors that can be constructed are then: 

1. Up{Pi)'Jf,Up{Pi) 

2. qyUp{Pi)a^^Up{Pi) 

3. {Piu + PAu)up{Pi)au^Up{Pi) 

4. [a] (Pi^ + PA^:)Up{Pi)up{Pi) 

N (P4M - Pl^.)Up{PA)Up{Pi) 

Out of these, (3) can be shown proportional to (4.b) by using Dirac equation, 
^Up{P) = iMpUp{P) as well as the usual 7-matrix algebra. It is also clear from 
eq. ( |2.9| ) that (1), (2) and (4. a) are linearly dependent. Hence (1), (2) and (4.b) 
exhaust the set of linearly independent fourvectors that can be constructed out of 
the spinors Up{P4), Up{Pi) and the fourmomenta Pi^P^ of the incoming and outgo- 
ing proton respectively. Hence the most general expression for the electromagnetic 
current of a spin 1/2 particle, which has a structure, can be written as. 



jr = '% 



n{q') Up{P,h,Up{P,)+n{q^) ^q^Up{P,)a,,Up{P, 



2Mp 



+Tl{q^) Up{P^)q,Up{P^ 



(2.10) 



where due to Lorentz invariance, JFf , jFf and ^f are three arbitrary functions of the 
scalars that can be constructed out of the available linearly independent fourvectors. 
Since the incoming and outgoing proton are on mass-shell, this means that Tl 
[i = 1, 2, 3) can be functions of g^ alone. Eq. ( f^.l(J| ) means that the structure of the 
proton can be parametrised by three arbitrary functions. 

We know, however, that the electromagnetic current is a conserved quantity, i.e., 
0. This translates in momentum space into 



d,JP^'^^' 



g^JP'-t = where q^ = P4>, - Pi^. 
Applying this condition to ( P^.1U| ) we get. 



f2.111 



= q,J'/ 



iqp 



Kr, 



J^{q )Up{P4)-f^qf,Up{Pi) + J^{q )^rT;rqp,qyUp{P4^)ay^,Up{Pi) 



2Mr, 



+ n{q')ur>{P,)q'up{P,) 



(2.12) 



The first term in the square bracket vanishes identically as can be seen by using Dirac 
equation for proton. The second term is identically zero since q^q^ is a symmetric 
tensor in (yU, v) and it is contracted with a^^ which is antisymmetric under (/i <-^ z/). 
Hence eq. (|2.12|) becomes 



n{q'Wup{P,)up{P^) = 



(2.13) 



This then tells us that Tl{q^) = 0. Hence the most general expression for the 
electromagnetic current of a spin 1/2 proton with structure, which is consistent 
with Lorentz and gauge invariance is given by 



Jr' = iqp 



Tl{q )Up(P4)7^Mp(Pi) + J^Kg )-^qyUp{Pi)a^^Up{Pi 



2M„ 



= Up{Pa)T ^Up{Pi) . 



(2.14) 



Thus we see that the most general expression for the electromagnetic current of 
an extended, spin 1/2 object is specified completely in terms of just two arbitrary 



functions J^i{q^) and J-'2{q^ 



I.e.. 



the effect of the spatial extension of the scattering 



centre on the scattering amplitude is simply parametrised in terms of these two 
functions. Comparing eq. (|2.14|) with eq. ( p. 81) we see that we recover the case for 
a proton with an anomalous magnetic moment Up but without a structure if jFf(g^) 
and J-^{q^) both are unity. Our discussions on nuclear form factors tell us that for 
q^ ^ l/(i?^)p, any form factor associated with the proton will be close to unity 



(recall eq. (|1.38| )). In particular eq. ( p.8| ) gives us. 



^f(O) 



Tm 



(2.15) 



The two functions J-'f{q'^) are called the two form factors of the proton and are 
associated with two linearly independent fourvectors in terms of which the electro- 
magnetic current can be decomposed. For a neutron which has no electric charge 
and only an anomalous magnetic moment we will have, 

J^{0) = 0; ^2"(0) = 1, (2.16) 

and Kp in eq. (|2.14|) will be replaced by k^. 



2.4 Electric and Magnetic Form Factors 

Having established the most general expression for the electromagnetic current of 
an extended object with spin 1/2, it is now necessary to calculate da/dVt for the 
scattering process of eq. ( |2.1| ) in terms of these form factors. We can calculate 
da/dQ by using eq. ( p.4|) for da, with A4 given by eq. (p.3|) , where J^ is given by 
eq. ( |2.14| ) and J^ by eq. ^]9|with a replacement of qp ^ qe, Mp — ;> me, P4 — ;> P3 and 
Pi -^ P2. We then get (see problems 4, Appendix B) 



%^''P ^ ''P^ 



p,ext 



-'mJ^J 11'-"^ -#-'': 



+^(-^f + «:p^fftan2^^. (2-17) 



In the limit q^ ^ 1/ {R^)p and with Kp = this reduces to eq. ( |2.5| ) as it should. 
Since the term (g^/2Mp ) tan^ | in eq. (2^) is known to arise from the magnetic 



moment of the "Dirac" proton and electron, on comparing eqs. ( |2.5| ) and ( [^.IT] ) we 
see that while the form factors jFf and J-^ were the natural ones while considering 
the tensor decomposition of JJl'^^^, the combination (jFf +fi;pjF|') has the more natural 
interpretation as the magnetic form factor. If we define the electric and magnetic 
form factors of the proton by 



Gl,iq')^J^!iq') + KpJ^^iq 



'\ 



GW)^mQ')-'^Pj[H^'2{q'). (2.18) 

p 



the expression for the differential cross-section becomes. 



da _ ( da\ p 

— - e p^ e p)= [-—] — 
dn V^^/MottPo 



+ TTTT^GIj {q ) tan - 



1 + gV4M2 2M2 ''' '^ ' 2 

(2.i9) 

The physical significance of G^(g^) and G\i[(q^) can be understood in terms of 
electromagnetic scattering of protons with definite helicity in the Breit frame. The 
latter choice means that the scattering involves no change in energy but only in the 



sign of the three momentum. To see this let us analyse the expression for Jp-*^^* of 
( ^.14| ), but for proton states of definite helicity (spin projection along direction of 
motion). Using eq. (|2.9|) we can rewrite this as 



^S = ^(A,A')r,n(Pi,A) 



iqpu{PA, A') 



.p,o, , ^{P., + Pl,){GlM')-Gl{q')) 



I.GIM') + 



2Mp (1 + g74M2; 



m(P4,A).(2.20) 



The Dirac spinor for a particle with helicity A when its momentum is along the z 
axis "pi = z\pi\ is given by 

u{puX) = nI \P1W3 I , (2.21) 

where 



X+1/2 and X-1/2 correspond to states with helicity +1 and -1 respectively. The Dirac 
spinor for a particle with helicity A, but moving along the negative z direction will 
be obtained by considering the transformation of the above spinor under rotation 
through 180° about y axis. This transformation is given by exp(— ivrcTj^), i.e., —iay. 
(The negative sign comes from the fact that we are rotating the physical state and 
not the co-ordinates.) Hence we get, 

/ -i(TyX\ \ 
ui-pi,X) = N \pi\a^ . (2.23) 

V^ + m/V 

In the Breit frame p^ = —pi and therefore we can use the above expressions for the 
spinors for a fixed helicity. We then get, 

^^;riA+i/2 = Gi.(a^^.(i^^^o,o), 

JTiI2,-iI2 = %GWW.-^)- (2.24) 



The first amplitude in eq. (|2.24| ) corresponds to no helicity flip (since the three 



momenta of the scattered and incident electron are opposite to each other, this 
means spin flip) scattering and this is proportional to the magnetic form factor, 
whereas the second amplitude represents the contribution with helicity flip (hence no 
spin flip) which corresponds to the electric form factor. Also note that the spin-flip 
amplitude due to the magnetic scattering vanishes in the non-relativistic (NR) limit. 
This observation also explains clearly the terminology used, viz., the nomenclature 
of electric/magnetic form factors. Thus electrostatic scattering cannot flip the spin 
of the electron whereas magnetic scattering does in this kinematical configuration 
(see discussion in problems). 



The discussion for 4§(e p ^ e p) will be applicable to the case of neutrons 



n ^>- e n] 



equally well. An expression similar to eq. ( ^.191 ) can be written for ^(e 

as well. Only the boundary conditions on GM,E{<f = 0) will change appropriately. 

Eqs. ( P^r^D , ( PTBD and (PT?1 ) imply, 



^^(0) = 1, GIM = 1 + Kp = fip (in n.m.) , 

GliO) = 0, Gl,{0) = Kn = fin (in n.m.) . (2.25) 

Experimental information on the form factors G^j^^{q'^) was obtained by studying 

^{ed -^ ed) and ^{ep -^ ep) = f ^j . The differential cross-section ^(en -^ en) 

is then obtained from the first two measurements applying correction factor for nu- 
clear physics effects. It is clearly not possible to do experiments with neutron targets 
as the neutron is not stable. Since we are discussing elastic scattering, eq. ( |1.23D 



ensures that there is only one independent variable characterising the final state. 



2 //I 2 ■ 2 /i/o\ / / 1 I ^^0 • 2 ^ 



This can be chosen to be either cos^ or g = (4pQsin 6/2) / 1 H sin 

Hence eq. (|2.19|) implies that 




( 



dnJep _ P ( ,p,^2, , DP^.2^^^^2^) ^ (2.26) 







where 



,2 



^n.^) - (TTWiMi) [^'^'(^'^^m-f^'(^'^ 




B'{<f) = ^ Gl,\q') 



Q 



2 



2M2 



(J-f + «:,J-f)^. (2.27) 



The above formula is called Rosenbluth formula. This also shows clearly how one 
can determine the two functions G^(g^) and G%{q'^) from a measurement of (j^j ■ 



From the kinematic considerations of eq. (|1.23|) , one sees that it is possible to obtain 
data at a fixed g^ and fixed angle 6, by changing the incident electron beam energy. 
The form factors for both the proton and neutron have been measured now over 
a wide range of g^ (g^ ^ 30 GeV^ for the proton and g^ ^ 5 GeV^ for the neutron). 
The experimentally measured form-factors for proton and neutron seem to obey a 
scaling law in the following sense: 



with Ml = (0.84)2 Q^y2 rpj^jg jg ^i^g ^jpQig fo^j^ f^p^Q^ of rp^^ig Hj^ ^^(g^) is 
identically zero as is to be expected. If we again compare ( |2.19D and (|2.5|), we see 
that to leading order in q'^, 

" -gV)^ .. . !,.... 4 • (2.29) 



( 



f )_„_.„. (1 + gVM^r 



ep, point 



This means that the elastic cross-section falls off very sharply with increasing values 
of g2 for scattering off an extended object. This behaviour then tells us that the 
probability of elastic scattering falls off very sharply with increasing g^ (z. e., increas- 
ing energy of the incident e~ beam). Some of the SLAC data on G\j{q^) taken from 
T5|| are shown in fig. ^ 

The experimentally observed g^-dependence of the form factors , eq. ( p. 291 ), 



along with eq. (|1.38| ) implies that the r.m.s. radius of the electric charge/magnetic 
moment distribution in a proton or neutron is 



(^■^>"'^ - owhv - '■'' f- • i^-^") 



The dipole form factor of eq. ( 2.29 ) can also be interpreted (see Table P as an 



exponential charge/magnetic moment distribution given by 

p{R) = -^exp{~MvR) . (2.31) 

Svr 

While it is true that the electric and magnetic form factors provide a neat way to 
parametrise the effects of spatial extension of the target on a scattering process, 
attempts to calculate these form factors from first principles, in some model, did 
not add to our knowledge of strong interaction dynamics or that of the structure 
of proton, beyond the information provided by eq. ( |2.31D . Real progress in the 
information on proton/neutron structure as well as strong interaction dynamics was 
made by increasing the energy of the electron beams incident on the proton/neutron 
target and studying the scattering process, which will be discussed in the next 
section. 

2.5 Summary of the dependence of the cross— sections for 
elastic scatterring on the nature of the target 

In this section we saw how the differential cross-section for elastic scattering depends 
on the nature of the scatterer and how this can be parametrised in terms of one (two) 
form factors for a scalar (spin 1/2) target. This information can be summarised in 
a compact form as in Table |^. 

As we can see, in the limit of the point proton, the form factors reduce to 1 
and the Rosenbluth formula reduces to eq. ( p.5|) as mentioned in the line 4 of the 
table. This does not incorporate the anomalous magnetic moment either. The term 
proportional to tan^ 9/2 in both these equations is due to the magnetic scattering 



Table 2: Dependence of the elastic differential cross-section on the nature of the 
target and projectile (Y: effect included, N : Not included). 



Formula 


Projectile e 


Target Z e 


Spin 


Energy 


Spin 


Anom. jj, 


Size 


Recoil 

(Target mass) 


Rutherford 
eq. (1.11) 


N 


Non Rel. 

(NR) 


N 


N 


N 


N 


Mott 
eq. (y^) 


Y 


Rel. 

E/m — > CX3 


N 


N 


N 


N 


eq. (^) 
spinless 
point p 


Y 


57 


N 


N 


N 


Y 


eq. (m 

point 
spin 1/2 p 


Y 


)) 


Y 


N 


N 


Y 


Rosenbluth 
eq. (^.IV) 


Y 


;7 


Y 


Y 


Y 


Y 



and is present when the spin 1/2 nature of both the target and the projectile is taken 
into account. As we have seen in our discussions at lower energies, the electrostatic 
scattering dominates and hence the second term can be dropped. If we drop this 
term from eq. ( p. 51) we recover the expression for the cross-section for the scattering 
of a spin 1/2 electron from a spinless, pointlike particle. For a spinles non-pointlike 
particle this would be multiplied by an appropriate form factor which at low energies 
of the projectile will reduce to one. If the target has no spin and has infinite mass (as 
compared with the projectile energy) then the target recoil in the scattering process 
can be neglected and i? ~ E'o for the projectile. In this limit this formula reduces 
to the Mott scattering cross-section of eq. (|1.42| ). At lower projectile energies, 
the effects of the spin 1/2 nature of the projectile also become negligible and the 
relativistic Mott scattering cross-section reduces to the Rutherford scattering cross- 
section. Again if the target is not a pointlike charge but a charge distribution instead, 
then this formula gets multiplied by a form factor as we have seen before. Thus we see 
that as we go to lower and lower projectile energies, various factors which represent 
effects of the projectile spin (e.g. the factor cos^(6'/2) in the Mott scattering cross- 
section of the second line in the table) or the target spin, size and mass ( the term 
with tan^ 6/2 factor, the form factor or the recoil factor respectively) go either to 
one or zero, giving us in the end the simple Rutherford scattering cross-section in 
the complete non-relativistic (NR) limit. 



3 Deep Inelastic Electron-Nucleon Scattering 

As we saw in the last section, study of elastic scattering of an electron off a proton 
target gives information about the spatial charge/magnetic moment distribution for 
a proton. However, eq. (|2.29|) and fig. ^ tell us that the elastic cross-section falls 



through four orders of magnitude as q^ changes from 2 — i> 25 GeV^. These larger 
values of q^ are reached using higher energy electron beams. With increasing g^, 
quasi-elastic scattering with excitation of baryon resonances becomes possible (recall 
the case of nuclei in nuclear reactions), and at still higher energies (and hence at 
higher q^ values), the scattering is dominated by the inelastic process: the so called 
Deep Inelastic Scattering (DIS). 

3.1 Kinematics of inelastic scattering 

The kinematics of the inelastic e~ - p scattering is shown in fig. ^, which is almost 
the same as fig. ^, except that the hadronic final state is no longer a proton. The 
reaction now is 

e- +p^ e- +X . (3.1) 

X stands for the hadronic final state. The four momenta of the various particles 
involved in scattering are as shown in fig. |^. The invariant mass of the hadronic 
final state is now, 

Ml = -Pl = -\p'\^ + W'^ . (3.2) 

Energy momentum conservation of eq. ( |1.6D gives us 

q^ = {P2 - P^f = (P4 - Pif = |p T -{W- M^f = -Ml + Ml + 2MpZ/ , (3.3) 

where u = Eq — E is the energy transfer from the electron to proton in the laboratory 
frame. Using eq. ( |1.9|) for g^ we then get, 

n 

q^ = AEEosm^- = -Ml + Ml + 2Mp{Eo-E) . (3.4) 



As we can see, if M^- = M^, i.e., the scattering is elastic, then eq. ( p^ ) reduces to 

g2 = 2MpU = 2Mp{Eo - E) . (3.5) 



Using the expression for q^ given in eq. (|3.4|) , we can see that eq. (|3.5| ) leads to 
eq. ( |1.23| ) with Ma — > Mp. Eq. ( |3.5| ) means that q^ and v are not independent 



variables and the reaction is characterised by only one independent variable. 

For quasi-inelastic scattering, i.e., excitation of a baryon resonance (say A^*"*" in 
the reaction e~p -^ e~N*^ — > e~p7r''), Mx = M* . Eq. (|3.3|) then becomes. 



e_ 

2 



g2 = 4^EEo sin^ - = -M*^ + Ml + 2MpV . (3.6) 



Again there is only one independent variable except that the relation between q^ 
and z/ is different from that for elastic scattering (cf. eq. (|3.3|) ). 



If Mjr does not have a fixed value but clianges continuously, then both g^ and z/ 
are independent variables. Clearly the proton no longer remains intact. This region 
is called the continuum region. This completely inelastic scattering is characterised 
by two independent variables. From eq. (|1.23|) and the definition of u it is clear 



that both of these are completely specified once the energy of the scattered electron 
E and its angle 6 are measured; i.e., the kinematics of the event is independent of 
the precise details of the hadronic final state X. It then makes sense to think of 
a measurement of the cross-section where one sums over all possible hadronic final 
states X. Such a measurement is called an inclusive measurement. 

The kinematically allowed region in the q^ - u plane for the elastic, quasi-elastic 
and inclusive inelastic scattering is shown in fig. ^ For elastic scattering the allowed 
region is the straight line given by eq. ( p.5|) . For quasi-elastic scattering, again the 
allowed region is a straight line but now with an intercept (M*^ — Mp)/2Mp, as 
given by eq. ( p.6|) . The discussions of last section, particularly eqs. ( p.l9|) and 
( p.29| ) , make it quite clear that the probability of elastic scattering goes down with 
increasing g^ as more inelastic channels open up. When both g^ and u are large then 
the cross-section is dominated by continuum excitation. The allowed region in the 
q^ - V plane in this case is the entire region to the right of the straight line given by 
eq. dp). 

The kinematics can be described in terms of any of the pairs of variables: {E, 9) 
or (g^, v). Equivalently, one can also define two dimensionless variables 

X = -^ = ~^' ■ (3 7a) 

2M.pV 2Pi • g ' ^ ' ' 

y = W = —^r- ■ (3.7b) 

For elastic scattering g^ and z/ are related to each other via eq. ( |3.5| ). Hence the 
variables x and y become, 

X = 1 ; 



2Eo/Mp sin^ ^ 
1 + # sin^ 



y ~ 1 , 2Bn -..2 e ■ (3-^ 



Mp " 2 



For inclusive, inelastic scattering both g^ and u (or equivalently x and y) can vary 
independently. We notice from eqs. (p.4[), ( |3.7| ) that 

g^ = 2MpUx = 2MpEoxy = -M|- + M^ + 2MpU 

= AEEosm^e/2 . (3.9) 

Since Mx > Mp, the equations above trivially yield the kinematically allowed region 
for variables {x, y) as 



Using eq. (|3l9| ) we can also derive the allowed region in the (g^, u) plane and it is 
given by 

The early DIS experiments at SLAC used electron beams with energy £"0 = 20 
GeV which corresponds to q^ < 40 GeV^. The current experiments at Fermilab 
use v beams with energy 500 GeV and hence can reach q^ values upto 1000 GeV^, 
whereas the e~ - p collider HERA at DESY operating with an electron beam of 30 
GeV and proton beam of 800 GeV is capable of measuring DIS cross-sections upto 
g2 < 10^ GeV^. 

3.2 Inelastic cross-section and structure functions 



We saw in the first lecture that for elastic scattering of eq. ( p.l| ), the most general 
expression for the current J^'*^^* and hence that for the cross-section could be writ- 
ten in terms of two arbitrary functions J-^{q^) and J-^iq^)- In case of the inclusive, 
inelastic measurement described above, again the cross-section can be parametrised 
in terms of two arbitrary functions Wi{q^, v) and W2{q^-, v) which have to be deter- 
mined experimentally. First let us see how the number of these arbitrary functions, 
which parametrise the effect of the structure of a proton on the DIS process, can be 
restricted to two using gauge invariance, Lorentz invariance, parity invariance and 
time reversal invariance. 

Let us start by recalling that in case of elastic scattering we started off with 
M. ~ {J^JJD/q"^, where J^ was given by eq. (|2.9|) and J^ for a point Dirac proton 
was given by eq. ( |2.9D whereas it was given by eq. ( ^.14[ ) for a proton with structure. 
In the present case, we need to sum over all possible hadronic states. Since now 
there is no spinor available to describe the hadronic state, unlike the earlier case of 
elastic scattering, it is not possible to repeat the steps which led us to eq. ( p.l4[ ). 



Recall, however, the calculation leading to eq. ( p3D , which gave the cross-section 
for an electron scattering off a point Dirac proton. Then we had, 

lA^Ve^V^;. where L,^ = ±Y.W\ K^ = Wr2 T. JPI' ■ (3-12) 

'^'^ S,S' ^^ Sp,S^ 

S, S'{Sp, S') denote the spins of the initial and final state electron (proton). In the 
case of a non-point-like proton all we had to do was to replace J^ of eq. ( p.9|) by 
jp,ext q£ gq_ ( |2.14| ). Comparing figs. ^ and ^, we see that the lepton end in both 
the elastic and inelastic case is the same. Hence we can once again write |A^p in 
a form analogous to eq. ( ^.12| ). However, in this case we know nothing about the 
tensor involving hadronic variables. 

To get more insight into the structure of the hadronic tensor for the case of 
inclusive scattering, it might be instructive to look at a few steps that led us to eq. 
31). The differential cross-section for the elastic process was written as: 

1 1 d^p 1 d^p' 1 e* 



da^ 



2Mp 2Eo(27r)3 2E {27i)^ 2W g^ 



J2 "e(P, S')-f^Ueipo, S)Ueipo, S)-f^Ueip, S') 



3,S' 



Op, Op 



X 



(3.13) 



The first factor in tlie square bracket is L^^ and tlie second is H'^^ of eq. ( p. 121) . 
The above can be rewritten as, 



da^ 



1 1 



cPp <fp' 



4 ^^^v 



2Mp 2Eo (27r)3 2E {2tt)^2W q 

Sp,Sp 
(2n)'6\P, + q - P,) 



with 



( P4, Sp\J^\Pi, Sp) = Up{p, Sp)jf,Up{pi, Sp) 



(3.14) 



Here J^ represents the electromagnetic current operator from which certain con- 
stants have been removed. The above equation can be used to write the expression 
for the double differential cross-section for elastic scattering as, 



^-^ tl.H'^ 



where 



H, 



1 1 



dVtdE g4 ^^ 
d^p' 



-'liu-' 



(3.15) 



E 



^"^ AnMp 2 ^ ., 

'' Op, Op 



(2nr 6\P, + q - P,) 



(27r)3 
(-Pl 5 'S'p I J I P4, S ) [S P4 \J^\Pi, Sp) 



(3.16) 



For the case of inclusive scattering we can generalise the expression for H^u given 
by eq. ( p.l6| ) as: 



H, 



1 



..-2i:i:/ii^n^(2-r*'(A+,-i: 



Sp X 

{P,,Sp\Jl\X){X\l\P,,Sp) 



4 = 1 



(3.17) 



where \X) denotes a state containing n particles with four momenta £±,£2, ■ ■ ■ An 
and we sum over all such states. The differential cross-section da*"^^ for this case of 
inclusive scattering is then given by (using eqs. (|3.13|) to (|3.16|) ) 



da 



inel 



2M'p Wo (2^2P ^ "^ ^^ ^^ 



(3.18) 



All we know about Hn^ given by eq. (|3.17| ) are its symmetry properties. Gauge 
invariance requires d^J^ = Q. In momentum space, this translates into 

q^H^, = H^,q, = . (3.19) 

We begin by writing H^y as the most general tensor consistent with gauge invariance 
and parity invariance that can be constructed out of the linearly independent four- 
momenta available at the hadronic end; g^ and Pi^. Also note that L^i, is symmetric 
in yU <-i> z/. Hence only the symmetric part of H^y will be relevant. The most general 
form for H^^y can therefore be written as 



H^y = W^{q\v)5^y + W2{q\ 



y 



p 



+W,iq^ uf-^ + W,{q\ u) ( ^^"'^^/"^^'- ) , (3.20) 

where Wi{q'^, v) (i = 1, 2, 4, 5) are arbitrary functions of q^ and v. The requirement 
of gauge invariance of eq. (|3.19| ) gives 

-.2 



<iv 



p p 



+ Plu 



q^ 



W^5(g^-)^ + 



M2 

p 



Pi-q 

M2 

p 



W2{q\v) 







Since q^ and Pi^ are two linearly independent fourvectors and the terms in the square 
brackets multiplying them are Lorentz scalars, it follows that each of the brackets 
must be identically zero. Hence we get, 

Pi-q. 



W,{q\v) 



q^ 



-W^iq^u) 



p p 



W^{q\v) = Q . 



(3.21) 



Using this, the expression for H^^ becomes, 
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(3.22) 

Thus we see that in effect the most general expression for the hadronic tensor involves 
only two arbitrary functions. L^j, of eq. ( p.l8|) is the same as that for the elastic 
scattering and is given by 
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-TrW^^Ifolu] 
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(3.23) 



To calculate da^'^'^^ of eq. ( |3.18D we need to know L^^H^y. Using eqs. ( |3.22| ) and 
(g), we get 

-4P ■ Po Wi{q\ v) + ^^^%f^ (2M2p ■ Po + 4P • PiPi ■ Pq) . (3.24) 
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With our choice of normahsation for Wj(g^, z/), (cf. eqs. (|3.18| ) and ( p.20| ), we have 



The superscript 'em' denotes here the nature of interaction involved in the scattering 
process, viz. electromagnetic interaction. In the laboratory Pi ■ Pq = —MpEo, 
P-Pi = -Mp{Eo - I/), and P ■ Po = -gV2 = -2EEo sin^ 9/2. Using these relations 
and eq. ( p.24| ), we get 






W2{q\u) + 2Wi{q\u)t^n^^- 



(3.26) 



Thus the cross-section for the inelastic scattering process e^p —^ e^X is parametrised 
in terms of two functions Wilq"^, v) and W2{q'^, v)- Since these functions contain all 
the information about the proton structure as revealed to an electromagnetic probe, 
these are called the structure functions. The differential cross-section of eq. ( p.26| ) 
can be equivalently written in terms of the pair of variables (g^, v) or the pair of 
dimensionless variables (x, y) introduced earlier. We can show that 

2MEox-—— = 2MEly—— = 2Mpn^ y- 



dxdy dxdq^ dq^di> E dVtdE 
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Before we go on to discuss the subject of experimental measurements of functions 
^i{q^i^)i ^2{q^,i^) and the information they yield about the proton structure, it 
is instructive to see to what the general functions Wi, W2 reduce to for the special 
cases of elastic scattering of an electron off a pointlike "Dirac" proton (eq. (|2.5| )) 
and a proton with structure (eq. ( |2.19| )). The cross-section of eq. ( |2.5|) can be 
rewritten as a double differential cross-section by using the identity, 

da"^ /■._./_ En \ da^^ 



,^ - , dE5\E „p" ,, — - . (3.28) 

dVl J \ 1 + i§ sin^ I j dfi ^ ' 

The above equation follows from realising that for elastic scattering the energy of 
the scattered electron is fixed via eq. ( |1.23[ ), once the angle is fixed. Using the 



definitions of v and g^, as well as properties of (5-function, we can rewrite the above 
expression for da'^^/dQ as 

'dE6 -u + ^] l + ^/sm^o -J?V ■ (3-29) 



dn J \ 2MpJ \ Mp 2) dn 

From this equation it is obvious that the double differential cross-section for the 
elastic case is 

-dndE= H~^ ^J V^k"" V ~dn ■ ^'-''^ 



Formally, eq. (|3.30| ) can also be derived by using 
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This is merely a restatement of the relation between E, Eq and sin^ 6/2 given by eq. 



Using eqs. ( |3.30| ) and (2^) we can therefore write, 
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Using the expression for {da/dQ)^^^^^ we therefore get 
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Similarly, we can see from eq. ( p.6|) that for the imaginary case of a spinless, pointlike 
proton we will get. 
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Using eq. (|2.19|) we get similarly. 
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(3.35) 



where A^i^q^) and B'^i^q^) are given by eq. ( [^.271 ). Comparing eqs. ( |3.33D - ( ^.35] ) 
with eq. (|3.26|) we get. 
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while 
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and 
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Table 3: The structure functions MpWi{x, Q^) and i>W2{x, Q^) for elastic scattering 
from different type of scatters. 



Scatterer 


MpW,{x,Q') 


i^W2{x,Q^) 


pointlike spin 1/2 proton 


4i;^(i-^) 


6(1 -x) 


pointlike scalar proton 





6(1 -x) 


spin 1/2 proton with structure 


j^^Gl,y5{l-x) 


/ \ 

1, -.V 



The functions G^(q'^) and G\j{q^) have the very steep q^ dependence given by eq. 
( p.29| ). This means that the structure functions W^f 2(0'^' ^) ^^^ ^^ very steeply with 
increasing q^. More interesting is the observation that the structure functions Wi^™" 
for a point scatterer do not depend upon the variables q^ and z/ separately, but are 
functions only of the combination x = q^ /2MpV. This is of course obvious. In the 
case of a point scatterer, the structure functions should depend only on dimensionless 
variables as there is no intrinsic length scale associated with the scatterer. Thus if 
the scatterer has a finite size then the structure functions for the elastic scattering 
Wf, W2^ , fall off as a power of g^ whereas for a pointlike scatterer they depend 
on g^ only through the combination q^ /2MpV. Table ^ summarises the behaviour of 
the structure function for different types of scatterers. 



3.3 Scaling of structure functions and partons 

The discussions of the earlier section tell us that the q^ and v dependence of the struc- 
ture functions Wi^2{q^i i^) depends on the nature of scatterer. Hence it is worthwhile 
asking, what do the experimental measurements of W^l '^{q"^, ^) look like? With in- 
creasing g^, the inelastic scattering begins to dominate the elastic process for g^ ^ 2 
- 3 GeV^ (recall eq. ( |2.29| )). For quasi-inelastic scattering, the resonance excitation 
corresponds to Mx = M* in eq. ( |3.4|) . Hence u is fixed once g^ is. This means that 
d?a/dq^dv will be a (5-function in v for a fixed g^. Experimental measurements of 
(Pa/dq^dv do indeed show peaks in i/ at a fixed g^. Fig. ^ taken from |]I6[ shows 
this. With increasing g^, the proton-resonance transition form factors also show a 
power law fall off with g^, just like G\i{q'^). At still higher values of g^, continuum 
production takes over. The interesting observation of Bjorken scaling was the fact 



that for this inelastic scattering the structure functions Wl^2^{q'^, v) do not fall off 
with increasing q^ but they are found to become independent of q^ (for a fixed value 
of q^ /2Mpv) as both q^ and u become large, i.e., they are functions of a: = q^ /2MpV 
alone. This is illustrated for some data on i'W^'^^{q'^., v) in fig. ^ taken from the first 
of Ref. B]. Thus the experimental observation is, 
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(3.39) 



This phenomenon of 'scaling' of structure functions was interpreted by Bjorken p| 
as an indication of the existence of pointlike scatterers inside the proton. This 
interpretation can be understood by recalling eq. ( |3.36| ). Thus the observed scaling 
of the DIS cross-sections indicates that inelastic electron-proton scattering can be 
understood in terms of incoherent, elastic scattering of electron off the individual, 
pointlike constituents of the proton termed 'partons' [Q. The charged partons which 
take part in the electromagnetic scattering are termed quarks. Whether these quarks 
are to be identified with the quarks whose existence is inferred from spectroscopic 
studies is best discussed later on in the context of the parton model. Since there is no 



scale associated with these pointlike objects, eq. ( p.36|) indicates that the individual 
elastic cross-sections must scale. We will show later, in a detailed discussion of 
the parton model, that the variable x can be interpreted as the fraction of proton 
momentum that the parton carries. Hence the intuitive picture is as shown in fig. 

|Tq|. 

In the scaling hmit of eq. ( p.39| ), using eq. (|3.27|) , eq. ( p.26|) becomes 
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S is the square of the total cm. energy of the scattering process ~ 2MpEo. The 
scaling limit corresponds to high energies of the incident electron, hence particle 
masses can be neglected and we get. 
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Here Fi^{x) and F2^{x) are the electromagnetic structure functions of the proton. 



It is worth noting here that according to eq. (|3.37| ), p^^^^caiar _ g Hence if 
the pointlike constituents inside the proton are scalars, the term proportional to 
tan^(6'/2) in eq. ( |3.26| ) (or equivalently the term proportional to xy^ in eq. ( |3.4CI| )) 
will be absent. Note also that all the discussions will be completely unchanged 
if one were to use /i~ beams instead of e~ beams. The structure functions are 
characterised by the target and the type of interactions used as a probe, in the 
present case a proton. So we have 
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3.4 Neutrino Deep Inelastic Scattering 

So far, in our discussions the scattering process that we considered was e~-p scat- 
tering. This probes the structure of proton via electromagnetic interactions and, as 
shown before, is characterised by two independent, arbitrary functions VFj(g^, z/) {i = 
1, 2) which are to be determined experimentally. One can also probe the structure 
of the proton via weak interactions in the charged current reaction 

ue+p^r+X (3.41) 

and the neutral current reaction, 

i^i+p^Ui + X . (3.42) 

Here we have started the discussion directly with the inclusive, inelastic process. The 
analogues of the form factors G^,^, G^^ and G^,j for the elastic processes involving 
neutrinos exist. Symmetries of strong interactions relate the electromagnetic form- 
factors as measured in the e~-p or /i~-p scattering and the weak form-factors as 
measured in i>£ p scattering. The charged current weak form-factors have been 
measured and indeed played a very important role in confirming these features of 
strong interaction symmetries. Here, however, we concentrate only on the DIS 



process shown in fig. |TI. 



Analogous to the earlier discussions (cf. eqs. ( p.l3| ) - ( |3.19| )) for the charged 
current reaction of eq. (|3.41|) , we can write 



^^"^'" = ..' Jlr. f^V (z^^X^-M,L'H'' (3.43) 
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where L' is the leptonic tensor evaluated using the weak current JJ['^'^'^{l,Vi) in- 
stead of the electromagnetic current used in eq. (|3.12|) , H'^^ is the hadronic tensor 
(again analogue of eq. (|3.17| ) but replacing the electromagnetic current by the weak 
current), g is the weak coupling of £, Vi to the weak gauge boson W and Mw is the 
mass of the W boson. The matrix element of the weak current is given by [^ 



(Ps, ^£| JriP2, Snu,) ^ Jr"(^, y^) ^ -^HSi, Pz)l,{l + l,)u.AS.a P2) ■ (3.44) 
Hence 
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Here, as before, masses of the leptons are neglected. We note now that L'^^ is no 
longer symmetric under an exchange /x <-^ z/. Hence the most general expression 
for H" must involve antisymmetric tensors too. A tensor decomposition of H" 
now involves six linearly independent tensors that one can construct out of the 
four momenta Pi and q, and hence six arbitrary functions of g^ and u. This is to 
be contrasted with the electromagnetic case where one needed only two arbitrary 
functions due to the its symmetric nature and the requirement of gauge invariance. 
The most general expression for H'' in this case can be written as 
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(3.45) 



Thus, in general, the DIS process z/£ +p -^ (i~ + X requires six structure functions. 
However, it should be noted at this point that the cross-section involves contraction 
of H'^j^ with L'^^, and contributions from W'^, Wl^ and W^ to the cross-section can be 
seen to be proportional to the lepton masses and can therefore be dropped. This is 
demonstrated for the term containing W'^{q^, v). Consider 
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Using -P3 ■ q = P2- q = q^/2] P^ ■ P^ = -g^/2 - ml we get. 
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Similarly terms proportional to W^ and W'^ in H''L' can be shown to be small. 

The contraction of the first two symmetric terms in eq. ( p.45| ) with the symmetric 
terms in V gives results similar to the electromagnetic case. The contraction 
between the symmetric and antisymmetric terms will obviously yield zero. The 
contraction of the antisymmetric term in H'' with the one in V gives. 
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2 [5aa'S(3f3' - Sal3'5f3a'] , WC get 



Using Pi ■ P3 = —MpE, Pi- P2 = —MpEf) and the expressions for P2- q, P^- q as well 
as g^ quoted earlier, we get 
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Neglecting the lepton masses, in the limit q^ <^ M^, the general expression for the 
inelastic, inclusive, differential cross-section for z/ reactions becomes. 
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where Gp/y^ = g"^ /SM'y^. The maximum i/-beam energy that has been reached 

in current experiments is Eq = 500 GeV. From eq. (|3.11|) this means that the 

maximum g^ that can be reached in these experiments is ~ 1000 GeV^. Since 

Mw^O (1000 GeV) , 

attainable values of g^. 

values of g^ are accessible only in the study of the reaction 



^ is valid even at these highest 
The charged current weak structure functions at higher 



the approximation g^ <^ M^r 
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u + X 



at the HERA collider where g^ can reach upto 10^ GeV^. 

Again, in the Bjorken limit of large g^ and z/ the structure functions are found 
to scale just as in eq. ( p.39| ), i.e., MpW^q"^ 
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S in the above equation stands for the square of the centre of mass energy given 
by Mp2 + 2MpEo ~ 2MpEo. Note that the result of eq. (|3^ ) is obtained from eq. 
( p.46| ) in the limit of vanishing particle masses just as in the case of eq. ( |3.40|) . Here 
we have three structure functions as opposed to the electromagnetic case. The third 
structure function arises from the term containing e^j^Q,^ in the tensor L'^^ {H'/^). 
This is the parity violating structure function. If instead of neutrino scattering we 
were to consider 



vi+p 
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+ X 



the corresponding term in L'^^, will change sign. Hence in the cross-section the term 
containing F3 will change sign. Hence, 



F^/ 



-E^- 



This can be physically understood by realising that z/(z/) is right (left) handed and 
the term containing F3 essentially arises from the V — A interference term in L' H'', 
which changes sign as we go from a right handed i^ to a left handed u.fj^ 

At this stage of analysis, it is not at all clear that the functions W[{q'^,h'), 
^2(5^5 ^) ^^d W^3(^^; ^) {or equivalently F^^{x), 2 = 1,3) have anything to do with 
the electromagnetic structure functions Wi{q'^,h'), i = 1,2 (or equivalently F^ 
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or Fl^^{x), i = 1,2). One can make predictions for relations between these only in 
the framework of the parton model. As a matter of fact, experimental test of these 
relations was an important step in establishing the parton model firmly. This will 
become clearer as we go on to discuss DIS and the parton model in the next section. 

3.5 Relationship between scaling of cross-section and point- 
like constituents 

The phenomenon of scaling of the cross-sections reflecting the existence of pointlike 
scattering centers occurs at different energy scales twice as we go from e~A to e~p 
elastic scattering to e~p deep inelastic scattering. This is illustrated by the data on 
nuclear scattering very nicely. Fig. O taken from ref. |]T3 shows this schematically 



whereas fig. |T3| fT^ shows some of the actual data. The "large" ness or "small" ness 
of a particular g^ value has always to be understood with respect to the inverse size. 
Recall, eg., that for q^ < (0.71) GeV^ the Wf^^ of eq. (g) will look just like their 



counterpart for elastic scattering from a pointlike scatterer of eq. ( p.36| ). In general 
therefore, for elastic scattering we can write, 

{M,Wf) uWf ^ f^,^2{x) g{q') . 

For g^ values such that g-Rtargct *^ 1 (^-c, q^ ^ ^target' where Atargct = inverse size 
of the target ~ (l/i?)), i'W2^ (MpWf^) will scale, i.e., will not show any extra g^ 
dependence. At these values of g^, g^q"^) ~ 1. In the case of elastic scattering, be it 
from a nucleon or a nucleus, the variables g^ and u always satisfy the relation 

g^ = 2Mtargot l^ ■ 

This also indicates that the function fi{2){x) must also have a factor 5 (jM^ — u 
or equivalently a factor = 6{x — 1), i.e., fi{2){x) will thus have a peak at x = 
g^/(2Mtarget J^) = 1- As g^ values increase, the function g{q^) starts differing from 1 
and falls off with increasing g^. For g^i?^ = 1, i.e., q^ ~ A^3^j.gg^, form factors cause 
a measurable suppression of the cross-section. For qR ^ 1, the function g{q'^) falls 
off very steeply indeed and the elastic peak at x = 1 disappears. 

In fig. |T|(a) the elastic peak at x = q^/2M^v = i at g^ = 0.01 GeV^ in eC -^ eC 
is shown. The quasi elastic excitation of the resonance C* in eC -^ eC* appears as 
a peak at a lower x values (cf. eq. (|3.3|)). At this value of g^, we have g^ ^ A^ and 
the Carbon nucleus appears like a point particle and elastic scattering dominates. 
As g^ increases further, at g^ = (0.1) GeV^, we have the situation 

A' oton > g' > A^ 



^proton ^=^ y -^ '^^C 



^^The terminology of left (right) handedness of the v(y) corresponds to its helicity being -1(1). 



In this regime the proton appears to be pointhke and the elastic peak for eC -^ eC 
disappears. The dominant process is no fonger the elastic process eC -^ eC, but the 
inelastic one which is now an incoherent sum of scattering off the N nucleons inside 
the Carbon nucleus. In principle the elastic scattering off a nucleon should show up 
as a peak at x = q^ /2Mc v = l/N (for the elastic ep scattering g^ = 2MpZ/; hence 
X = q^ /2Mc V = 1/iV)- The Fermi motion of the proton in the nucleus however 
changes the kinematics and smears the (5-function peak. This is shown in fig. 0(b) . 
Now 

But g'^{q^) ~ 1 since q^ <^ ^proton- Thus /^(x, Q^) now scales. This scaling thus 



reveals the existence of A^ pointlike nucleons in the nucleus C. Fig. |13] shows data 
on ea scattering at two different q^ values q^ ~ 0.08 GeV^ and q^ ~ 0.1 GeV^. This 
figure illustrates the same point, only the value of N here is 4. 

If we now consider a proton at rest, at higher values of g^ (~ 0.5 GeV ^) and 
higher v values, the elastic peak for ep — > ep will occur at a; = -^ and the quasi- 
elastic excitations corresponding to ,say, A++ will show up as a peak at still lower 
values of x. This is the proton analogue of Fig. |T2|(a) and which was discussed 



already in section 2; except for the fact that the variable defined there had been 
in terms of proton mass {xp = (g^/2MpZ/)). Now the cross-section (more precisely 
i/Vrf) again will consist of two factors: one an x dependent kinematic function and 
a q^ dependent form factor reflecting the target size. This is the violation of scaling. 

Now if g^ is further increased to g^ ^ 0.71 GeV^ (Fig- |12|(d) corresponds to 
g^ ^ 5 GeV^), then the rapid fall-off of the form factor, g'^iq"^), with increasing g^ 
will cause the elastic peak to vanish as well as the inelastic A"*""*" peak. Now the 
incoherent scattering from the pointlike partons (quarks) eg — > eg will show up as 
a smeared peak in x distribution. The x value at which this peak appears should 
give information about the number of partons off which the electrons get scattered. 
The structure function scales again, revealing the existence of pointlike constituents 
inside the proton and this peak will be at x ~ 1/3N indicating the existence of three 
valence quarks. 

At still higher g^ values, for the proton target, (g^ ~ 200 - 400 GeV^ for the ep 
, fip and up experiments discussed earlier or g^ ^ 10^ - 10^ GeV^ at the ep collider, 
HERA) gg production begins and that increases the number of pointlike constituents 
in the proton effectively. This shifts the peak in the structure function to lower and 
lower X values. These scaling violations are better discussed in terms of QCD and 



have been studied extensively in current DIS experiments ||19|| . However, these will 
not be discussed here any further. 

4 Parton Model 

4.1 Formalism 

As discussed in the last chapter, the DIS e~p ^ e~ X cross-section scales at large 
g^ and large v (in the Bjorken limit). The observed scaling is evidence that the 



DIS cross-section is given by an incoherent sum of scattering of the electron against 
individual partons inside the proton. This interpretation is basically the parton 
model. Implicit in this model are two assumptions: 

• Interactions among the partons are negligible during the time of interaction 
between the electron and the parton. The higher the energy, shorter is the 
interval of time of this interaction. Before the proof of asymptotic freedom 
P of QCD, this assumption was justified only by the success of the parton 
model, but asymptotic freedom provides now justification for this assumption. 

• The second assumption is that final state interactions can be neglected. If the 
struck parton receives a huge kick then it gets removed from the parton and 
the final state interactions are hence negligible (recall here fig. p!0|). 



Fig. n shows a schematic description of DIS in the parton model. Let ^ be the 
momentum fraction of the proton carried by the struck parton and Cq denote its 
electric charge in units of the proton charge. Let fi/p{C,) be the probability that 
the parton i carries momentum fraction ^ of the proton. Momentum conservation 
implies 

Y. f dU■^/p{0 ^ = ^ ■ (4.1) 



Note here that the sum is over all types of partons, not just the charged ones which 
the incident electron (or equivalently the virtual photon) sees. In this model, one has 
neglected transverse momentum of the partons. This is justified by the experimental 
observation that apart from the struck parton which causes hadrons to emerge at 
large angles, the remaining particles in the final state emerge at small angles w.r.t. 
the beam direction. 

Let us choose the z axis to be the direction of motion of the proton and hence 
of the parton. The magnitude of the three momentum is \pi\ and C,\pi\ respectively 
for the proton and the parton. The four momenta are given by Pi = {0,PL,iEi) 
and Pq = (0,^Pi,zC^i), where Pl = \pi\. This gives Pf = -M^ and P^ = -^^M^. 
It appears as if the partons have a variable mass ^Mp. This clearly is not what we 
mean. This is just a reflection of the fact that the above kinematics is strictly correct 
only in the limit where all masses can be neglected. In this case the kinematics given 
above simply corresponds to a collinear emission of a massless particle from another 
massless particle. 

The frame of reference in which the above kinematics is strictly valid is called 
the infinite momentum frame. In the infinite momentum frame, time dilation slows 
down the rate at which partons interact with one another and this time scale is now 
much bigger than the time taken by the current {i.e., the electron or the virtual 
photon) to interact with the parton. Hence the impulse approximation (assumption 
(i) above) is justified in this frame. 

After scattering the struck parton (quark) appears in the detectors as a stream 
of hadrons. The time scale of hadronisation (~ 10~^^ sec.) is much bigger than 
the interaction time scale ( ^ 10~^^ sec.) for u ^ 0(10 GeV). This description of 
scattering as a two step process is the second basic tenet of parton model. It is clear 



from the above discussion that for both these assumptions to be justified and the 
picture of fig. |T1| to be true, g^, z/ and W all need to be large. 

Since the partons are pointlike objects, mgiyf'^'^*°°(g^, u) and 1/1^1"^'^*°° (Q'^j ^) ^.re 
given by eqs. (^.361), replacing Mp — >• nig = ^Mp. Hence we have. 
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If we now define 



then we get 
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The right hand side of this equation is clearly a function of x alone, and we can 
write (with an analogous discussion for vW2^^^°^{q^, z/)). 



Fr'°"{x) = iyWr'°^q^ ly) = 6 (l - j) 
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Since in the hmit of large g^, u and W, the total cross-section for ep scattering (and 
hence MplVf^"^ °°, uW2'^^^°^) are given by an incoherent addition over all the charged 
partons (quarks) we get in the scaling limit. 
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Similarly, 



frw = ii:/>.;^^« -.)/,/,(?) 



It should be noted here that the summation is over particles as well as antiparticles. 



Thus we find that the scahng variable x which we defined earher can be identified 
with the momentum fraction ^ that the parton q carries. The structure functions 
appearing in eq. ( |3.39| ) (which are the scahng limits of the arbitrary functions 
Wi{q'^,h'), i = 1,2 which appeared in the tensor decomposition of eq. ( |3.22D ) are 



related to the probability of finding a parton of charge Cg (in units of proton charge) 
with momentum fraction x of the proton. There is yet another way of seeing the 
same thing. In the infinite momentum frame all masses are negligible. The final 
four momentum of the struck parton is given by {C,Pi + q)- Hence we have 



which gives us, 



(^Pi + g)' ~ 



-r 



(4.4) 



2Pi-g 

This is same as the variable x defined in eq. (|3.7a|) . The above expressions for 
F^^{x), i = 1,2 also imply 

Ff (x) = 2xFf (x) , (4.5) 

which is known as the Callan-Gross relation. Using eq. ( [4. 5) ) we get from eq. ( |3.40D 
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(4.6) 



In the parton model picture which we have developed above, now the structure 
function F2^{x) is expressed in terms of the probability density functions fq/p{x). 
Since this probability should be independent of the probe that is used to extract it 
from the data, one expects that structure functions F^^{x) {i = 1,2,3), measured 
in up DIS, must also be some combinations of the same functions fq/p{x) that are 
extracted from the data on electromagnetic DIS processes, the specific form being 
decided by the nature of uq interactions. As a matter of fact, identifying the charged 
partons in the proton with the constituent quarks of SU{3) fiavour, one can make 
definite predictions for the weak structure functions of the parton as well as elec- 
tromagnetic and weak structure functions of other targets such as neutron or nuclei 
and relate them to each other. An experimental verification of these relations in the 
u DIS experiment |^ played a very important role in establishing the parton model 
on a firm footing. 

An alternative definition of eq. ([4.6|) can be given as follows. Consider the cross- 



section da/dVt for scattering of an e from a pointlike object of charge Cq (in units 



of proton charge). The corresponding expression is given by eq. ( |2.5| ). For the case 
of elastic scattering, we know that y is given by (recall eq. ( |3.8| )) , 
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Hence 
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where S is the square of the cm energy. In the present case we have to consider scat- 
tering of the e~ from a charged parton carrying a fraction x of the four-momentum 
Pi of the proton. The square of the cm energy S of the elastic electron-parton 
scattering is given by, 



S = -{xPi + Ps)^ = -2xPi ■ Pa = 2xMpEo = xS 



(4.7) 



Using the expression for da'^^ /dVt given by eq. ( |2.5| ) but now for cm. energy y S" as 
given by eq. (^^), we have (again neglecting particle masses), 
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If fq/p{x)dx is the probability of finding a parton with momentum fraction of proton 
between x and x + dx, we have 
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dy 

Hence we get 
da 
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Note that above equation is the same as eq. ( [4.6|) or eq. ( p.40|) where 
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(4.10) 



This alternative derivation of eq. (4^) indicates that, to derive the corresponding ex- 
pressions for the double differential DIS cross-section for vp or vp processes, we need 
to know da"''/dy{h'e q ^ i q'), da"''/dy{h'e q ^ i q') etc. Due to the parity violating 
nature of weak interactions, the angular distribution (and hence y distribution) of 
uq and uq scattering are different in nature. This is in contrast to the situation in 
the case of electromagnetic scattering. The differential distributions da'^^/dy for the 
elementary scattering process can be obtained after a simple calculation to be 
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The second equation above indicates the impossibihty of vi q (z/£ q) scattering in the 
backward direction {6 = it corresponds to y = 1). This can be easily understood 



from fig. |T^. This shows that the backward scattering in this case will correspond to 
lAJ^I = 2 which is not possible. We can use the above expressions for da^^ /dy{vi q) 
etc. in a manner similar to the one used in ariving at eq. ( ^.9| ). In this case we get 
for the charged current DIS cross-section, 
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Of course it is understood that the sum is to be taken over those quarks or antiquarks 
which can take part in allowed transitions. The first term in the bracket is parity 
conserving whereas the second one violates parity. This equation is the same as eq. 



( p. 47] ) if we identify 



F^'^ix) = 2xFl''{x) = 2x 



Y. fi/pi^) + fg/pi^) 



(4.11) 



xF. 



up I 



X] 



2x 



Y fg/pi^) - fg/p( 



I q 



X) 



(4.12) 



4.2 Structure functions for proton, neutron and isoscalar 
targets 

In this section we will write the form for the electromagnetic and weak structure 
functions expected in the parton model, if one identifies the charged parton with 
the Gell-Mann-Zweig constituent quarks. In the SU{3)f picture, e.g., the proton 
contains two w-quarks with charge |ep and one ci-quark with charge — |ep . Hence 
one expects that F|^ is given by. 



F,''' 
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-xu^{x) + -xd^{x) . 
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(4.13) 



Here u^{x) = fu/p{x) and so on. Isospin invariance would imply, 

m"(x) = d^ix); rf"(x) = u^ix) . 
Hence we expect. 
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U^ix) + ^u^ix) 



(4.14) 



However, this presupposes the picture that the proton (neutron) has 2(1) w-quarks 
and 1(2) ci-quarks. However, all the conclusions, verified by experiments, about 
static properties of proton/neutron will remain unchanged, if in addition to these 
partons one had large number of uu, dd, ss pairs forming an SU{3) singlet, which 
have been radiated from those "valence" (large momentum) partons. One would 
expect the proton/neutron to contain such radiated partons in the parton model 
picture. Since these arise from radiation, these will have carry less momentum than 
the parent quarks, they will have a softer momentum distribution than the "valence" 
quarks. These were called "wee" or "sea" partons by Feynman. Then we will have, 
assuming SU{3) symmetry for the sea {i.e., the sea quark density is the same for all 
three types of quarks u, d and s), 



U^{x) = Uy{x) + U^si^)y 

rfP(a;)=4(x) + 4(x), 
u^six) = u^si^) = dl{x) = (Fs{x) 
In this case eq. ( [4.13| ) will modify to 
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Again isospin invariance gives us 
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In writing eqs. ( ^. 16| ) and ( ^4.17| ) we have assumed, in addition to the fiavour sym- 
metry of the sea, also absence of heavier charm and bottom quarks in the sea. At 
higher energies, even these can be radiated. In that case the factor before K{x) will 
change. Using eqs. (|4.16|) and ( [4.17|) we get. 
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(4.18) 



Can we compute these structure functions u\r{x), d^{x) and K{x)l The answer 
is no, not without a model. However, we can try to measure them experimentally. 
The experimentally measured DIS cross-sections for ep (fip) reactions will yield 
F2^{F2^). How does one extract the specific parton densities from it? Before one 
goes into this question, let us see how one can obtain some qualitative information 
x) and K{x). Consider the ratio of eq. ( [4.18|) . If K{x) dominates 



about Uy{x 



dU 



over Uy[x 



dvix) 



this ratio will be 1. If either My(x) 



or dy{x) 



dominates over all 

the other densities then this ratio will go to 1/4 or 4 respectively. Hence, under the 
assumption of valence dominance we will get 



1 F^'^ix) 
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(4.19) 



If, on the other hand, we have K{x) = and Uy{x) = 2dy{x) (as can be expected 
in a pure constituent picture), this ratio will have the value 2/3. 

With these qualitative predictions we turn to the question as to what the ex- 
periments say. As was the case with the form factors, F2^{F2^ 
combining the data on ep, ed {fip, fj,d) scattering. Early data[^ tell us: 



2 ) is measured by 



fa) K(x) dominates at small x and the ratio is indeed close to 1 at small x. 



(b) At large x{x -^ 1), F2^{x)/F^ 



[X 



4. This means that Uy{x) dominates 



over dy{x), as well as over K{x) at large x. 

Observation 
not fuUfilled. 



(c) Observation (b) also tells us that the naive expectation of Uy{x) = 2dy{x) is 



Fig. ITB| shows the recent high statistics data of the NMC collaboration on the ratio 



p^n jp^p [^. What are the other theoretical constraints on these densities? We 
know that the net number of u{d) quarks in a proton is 2(1) and the net number of 
strange quarks in both is zero. Hence we have, 

{nF{x) — u^{x)) dx = 2 = / Uy{x)dx , 
Jo 

(d^(x) — (F{x)] dx = I = / dy(x)dx , 
[SP{x)-SP{x))dx = . 



(4.20) 



The first two constraints of eq. (f4.20|) can also be obtained by considering elec- 
tromagnetic charge conservation, which implies, for the proton and neutron respec- 
tively. 



1 
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dx 
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(4.21) 



We have to use the fact that both the proton and neutron have zero strangeness, 
in order to get the third of eq. ( |4.20| ). It should be emphasized here that the 
actual determination of uP{x), d^{x) and K{x) has to be done by combining data 
on electromagnetic structure functions as well as the weak structure functions for 
different targets; p, n and isoscalar nuclei and then fitting a form to the densities 
Uy{x), dy(x) and K{x), using the data, subject to the above sum rules. Some of 
these details will be discussed below. 



4.3 Properties of partons as determined from DIS 
4.3.1 Spin of partons 



Let us start from the Callan-Gross relation of eq. (15). Note here that if the charged 



partons were scalars then, according to eq. ( |3.37| ), Fi^'^"^(x) = 0. Hence the ratio 
F2^{x)/{2xFi^){x) -^ oo whereas eq. (|475| ) predicts the ratio to be 1 for spin 1/2 
partons. It is customary to define a longitudinal structure function F'^{x) by 

Ff (x) = Ff (x) - 2a;Ff (x) . (4.22) 

Since F2^(x) and Fi^{x) can be extracted from the measured DIS cross-section, an 
experimental verification of the above relation will imply that the pointlike con- 
stituents inside the proton revealed in DIS have spin 1/2. F^{x) is termed the 
longitudinal structure function as it is proportional (in the large g^, v limit) to the 
virtual photoabsorption cross-section for longitudinal photons {i.e., photons with 
helicity A = 0). In fig. |1^ we see the experimental data as the ratio 

Fr^ix) 
2xFi{x) 

In the large g^, z/ limit this is the ratio of virtual photoabsorption cross-section for 
longitudinal and transverse photons (A = ±1). As noticed before, the denominator 
is zero for scalar partons. This can be physically understood as the impossibility of 
absorption of transverse photons by a scalar target. The data of fig. [l^ taken from 
pl[| show clear evidence that the charged partons are spin 1/2 objects and not spin 
objects. Deviation of this quantity from zero is yet another 'check' for QCD but 
again will not be discussed here further. 

4.3.2 Momentum carried by charged partons 



As can be seen from eqs. ( |4.16|) and ( [4.17|) the area under the F2^{F2"') vs. x curve 



will measure the weighted sum of the momentum fractions carried by the charged 
partons in the proton (neutron). Hence 

/o F2^{x)dx = |e„ + |ed + |e, , 
^ (4.24) 

where e„,erf and e^ are the fractions of the proton momentum carried by {u + u), 
{d + d) and (s + s). If we define 5 = es/(e„ + e^) , we can write 

£ [Ff (x) + F^ix)] dx = ^(e„ + 6,) + p{e^ + e,) , 

^^^ C [i^r (^) + i^r(a^)] rfx = e. + e, + e, . (4.25) 
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It is clear from eq. ( [1.25| ) that a knowledge of 5 is necessary to determine (e^-l-es-l-erf) . 
An extraction of quark densities from the DIS data using ep, en, up, un experiments 
and eq. ( [4.28| ) shows that S < 0.06. Using the data on {F2''{x) +F2"'{x)), this gives, 

Cu + Cd + e, ~ (0.54 - 0.56) ± 0.04 . (4.26) 



Eq. ( |4.26|) above implies that some momentum of the proton (~ 50%) is carried by 
partons which are neutral to the probe, i.e., partons which do not have electromag- 
netic or weak interactions. These are called gluons which hold the quark-partons 
together in a proton. 



4.3.3 Charge assignment of different partons 



To confirm the fractional charge assignment used in arriving at, e.g., eqs. ( [4. 131 ), 
( [4.14| ) one has to combine data on F2^, Fg" with the data on F2^, F2"'. For this we 
should obtain expressions for ^2^^, Fg^" analogous to eqs. ( [4.13| ), ( |4.14| ). To do this 
one has to note the following : 



(i) A neutrino can scatter off only the charge — |ep quarks and charge — fcp anti- 
quarks, as it has to have a i~ in the final state. Hence only Ui d, vi s and vi u 
processes take place. (Recall, we have at present neglected the heavier, charm 
and bottom, quark content of the proton/neutron). 

ii) At low energies, Ui s ^ i'c transition can be neglected and z/^ s -^ £~u 
transition is suppressed in the limit of the mixing angle in the s-d sector 
(Cabibbo angle) 9^ — 0. This is a good approximation as cos^c = 0.98. Under 
these approximations, the only transitions that contribute to the DIS cross- 
section for incident z/^ are ui d (u) -^ i~u (d) and U£ u (d) -^ £^d (u) for 
incident ue. 



Using eq. ( |4.12|) then we have, at low energies and with the above approximations, 



F^^'ix) = 2x [d P{x) + mP(x)] 
F2^"(x) = 2x\uP{x) + dP{x) 



(4.27) 



Using eqs. (^IBj) , (^Tzj) and (|]23) we get 



F^Pjx) + F|"(x) 



uP{x) + d P{x) + uP{x) + dP{x) + SP{x) + SP{x) 



uP{x) + dP{x) + uP{x) + dP{x) 



(4.28) 



Excess of this ratio above 5/18 is a measure of the momentum carried by the s{s) 
quarks. Production of charmed particles in DIS in vp reactions provides a direct 
measurement of this strange sea content of the proton. The factor of 5/18 in the eq. 
( [4.28| ) above is the average squared charge of the m, d quarks. 

The electromagnetic and weak structure functions are usually measured not only 
for light targets such as proton/deuterium, but more often experiments are per- 
formed with heavier, nuclear targets so as to get large cross-sections. Normally 
one expects that for nuclear targets, the cross-section from different nuclei will add 
incoherently. For an isoscalar target therefore, the structure function per nucleus 
becomes. 
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With eqs. (|]T|) and ( pTj) , we get 
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If we include also the charm sea contribution then we will have 
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Note again the factor of 5/18. It is clear from this discussion that the ratio 

the same as the r.h.s. of eq. (|4.28| ). A collection of some of the data taken from 

the third of Ref. [^] presented in fig. |T8| shows experimental proof of this relation. 

This confirms the identification of the charged partons with the Gell-Mann-Zweig[jl| 

quarks. 

/•I 

Note also that a measurement of ^ / F2 dx for an isoscalar target can directly 

give information about the total momentum carried by the quark-partons. The 
recent high Q^ data on fip/^d DIS scattering by the different collaborations like 



EMC, BCDMS, NMC and up/uFe DIS data by the CCFR collaboration |T9[ confirm 
all the parton model features quite beautifully. 



4.4 Interpretation of sea densities 

The existence of gluons gives a very simple understanding of the sea-quark den- 
sities. At some very low g^-scale, the proton can be looked upon as made up of 
only three valence quarks u, u and d, all carrying (1)"^*^ of the proton momentum. 
The probability distributions Uy{x) and dy{x) are just 5{1 — x/S) with appropriate 
normalisation. Emission of gluons, by bremsstrahlung, by these quarks causes the 
probability distributions to shift to lower x values. The emitted gluons give rise to 
qq pairs thus giving rise to sea quarks. Since the process of bremsstrahlung is natu- 
rally peaked at small values of momentum fractions, it generates qq sea. At higher 
and higher energies the number of qq pairs produced goes on increasing. Actually 
this process will also give rise to scaling violations as it causes vP{x), d^{x) to shift 
to lower X values with increasing g^ and makes F2^{x) q^ dependent. However, this 
clearly takes us out of the realm of the Quark- Parton- Model (QPM) and causes cor- 
rections to the simple QPM picture. These corrections (scaling violations) actually 
played an important role in establishing the nature of interactions among quarks 
and gluons. But this will not be discussed further here. 

Another way of understanding more about the sea-quark densities is to try to 
construct a quantity which is independent of sea quark densities. Consider the 
following combination: 



X 



Fl'{x)-F^-{x) = -[uUx)-d'v{x)] 



(4.30) 



This difference does not involve the sea quark densities at all if one assumes isospin- 
symmetric sea densities. The experimental data on {F2^{x) — F2^"(x)) [BTI] show a 



peak at X = 1/3. Fig. |T^ shows (just to show the increased accuracy of the newer 
data) the much more recent and high statistics data taken from [^. This clearly 
supports the picture of three constituents of mass |Mp and also the interpretation 
of sea quark pairs as arising due to bremsstrahlung from valence quarks. 

4.5 Sum rules on parton densities 

Since there is, as yet, no theory which can compute parton densities in a proton, 
all the knowledge about parton densities is to be obtained from experiments, which 
can then be used for testing models of the proton structure functions. Within the 
framework of the quark-parton model, various relations, sum rules have been derived 
for Fg^ and F2"' or combinations thereof. Some of these sum rules have already been 
written down in eqs. ( [4.2(]| ) and (|4.21|) . From eq. ([4.30|) we get 



/' i^l^^JUMdx = \f K(x) - 4(x)) dx = \ . (4.31) 
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This is called the Gottfried sum rule. This sum rule is arrived at by assuming 
vPs{x) = d^{x). If that is not the case, this will be violated. Current data from 
H DIS experiments |2^ show that the sum rule is violated, and the above integral 
(obtained from extrapolation of the data for the region x > 0.8 and x < 0.004 ) is, 

1 F^'^ix) - F^^ix) 

^ ^ ' ^-^^t/x = 0.258 ± 0.017 . 

X 

These data use measurements of the structure function at g^ = 4 GeV^. These 
indicate a departure from the expected value of 1/3 and a breaking of the isospin 
symmetry of the sea densities which has been assumed in arriving at eq. ( [4.31| ). 

One can derive yet another sum rule using F2^ and F!^"'. We see from eq. ( |4.27|) 
that , 

"1 (^2^" - F!^^)dx n 



2 / {uP{x)-d'P{x))dx = 2 . (4.32) 
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This sum rule is called the Adler sum rule. Again it assumes isospin symmetry for 

the sea-quark densities. The sum rule is known to be satisfied. But the z/ scattering 

data have intrinsically much larger errors compared to the fip{^D) scattering data. 

As a result the neutrino data are not in a position to test the violation of isospin 

symmetry of sea densities implied by the /j, DIS data from EMC. 

Yet another sum rule that has been written down is for the parity violating 

structure function F^^. One can write down expressions for F^^ and F^"" using 

eq. ( ^.12| ) and using the isospin invariance of the neutron/proton parton densities. 



Using eq. (|4.12|) we get, 

j\Fr{x)+Fr{x))dx 



This is called the Gross-Llewellyn Smith sum rule. The earliest v experiments |]2^ 
verified this sum rule. Currently it has been a focus of lot of discussions as the 
deviation of the r.h.s from three can provide important information about and test 
of pQCD. The currently measured value is 2.50 ± 0.018(stat.) ± 0.078(s|/st.). 

4.6 Parton model in processes other than DIS 

DIS processes provided the first measurement of parton densities. We saw in the 
previous sections that a combination of the different (weak and electromagnetic) 
structure functions with different targets (proton, neutron, isoscalar nuclei) can be 
used to establish different properties of the charged partons such as their spin, charge 
assignment etc. The individual parton densities can be extracted from the data on 
F2 only under certain assumptions and one usually fits a form to these. The fits are 
constrained by various sum rules which are derived on general principles. However, 
none of the DIS processes can give information about the gluon densities. In the 
framework of perturbative QCD, outside the realm of quark-parton model, some 
information about gluon densities is obtained by studying the scaling violation (g^ 
dependence) of the structure functions. But no direct information on the gluons is 
available from these processes. It is also important to note that since there are, as 
yet, no theoretical predictions for either quark/gluon densities, it is imperative to 
find processes other than DIS to extract combination of parton densities different 
from those measured in DIS, to supplement our knowledge. This, in turn, can help 
to get a more complete picture of parton densities. 

A better knowledge of parton densities is not only essential to check our ideas 
about perturbative QCD but it is also necessary to be able to make correct pre- 
dictions for the cross-sections of different physical processes expected at the high 
energy e-p or hadron-hadron colliders. High energy collisions of hadrons (say p-p, 
TT-p or TT-A) with each other can be described in parton model as an incoherent 
sum of interactions among partons. This means that our predictions as to what is 
likely to happen at high energies in collisions of these hadrons will depend on our 
knowledge of these parton densities. 

As an illustration of the above, consider the production of a vector boson W^ 



in pp collisions as shown in fig. |20|. The production cross-section is given by 



da{pp -^ W^x) = / dxi / dx2nF {xi)<F {X2) da{ud — ^ W^) 

(4.33) 
The important point here is that the densities u^ or # in eq. ( ^.33] ) are exactly 
the same as those extracted from DIS. That the same probability functions are 
applicable in all 'hard' processes be it DIS or be it (say) W'^ production, is an as- 
sumption in the QPM. In the framework of perturbative QCD this has been proved 
quite rigorously for a large number of hard processes]^]. So any hard processes in 
hadron-hadron collision is computed as an incoherent sum over all the partons where 
the individual 2 ^ n parton subprocess is convoluted with the parton distribution 
functions. An agreement of measured cross-section with the predictions made in 



the QPM provides consistency checks on our knowledge of these parton densities. 
Further, by singhng out final states which are sensitive to a specific parton in the 
initial state, we can better our knowledge of parton densities. For example, a com- 
parison of W^ production discussed above with W~ , can be a good test of the 
isospin symmetry of the sea. 

This simple QPM picture has been justified theoretically in perturbative QCD, 
has been used extensively and has strong experimental support. The high energy 
processes of interest to the Nuclear physics community are the heavy ion collisions. 
The plain QCD contributions to different final states in these collisions, arising from 
the partons in the nuclei, need to be computed correctly before one can assess the 
observability of different signals of Quark-Gluon-Plasma (QGP) formation such as 
J/ip suppression, dilepton production or direct photon production [^. For this one 



needs a good knowledge of parton densities inside the nuclei. This is the topic of 
the next lecture. 

4.6.1 Jet Production 

The most common hard process in hadronic collisions is called jet production. This 
arises basically from the scattering of two partons in the hadrons which produce light 



partons in the final state. As already said in eq. ( [4.33| ), in QPM the total cross- 
section at the hadronic level is obtained by convoluting the parton densities with 
the cross-section of this basic 2 — »• 2 subprocess. Since this scattering is between 
pointlike particles, it involves large momentum transfers and the final state partons 
are thus produced at large angles {i.e., with large momenta transverse to the beam 
direction, {pr)) as opposed to the remaining partons in the two hadrons which go 
along the beam direction. There are in all 8 hard scattering 2 — i> 2 subprocesses 
which will produce light partons in the final state. The g or ^f in the final state 
does not appear in the detector as a free gluon or quark due to the phenomenon of 
'colour confinement'. They appear in the detector as a 'stream' or 'jet' of particles 
clustered around the original direction of the parton. As said earlier, the major 
feature of QPM is that these jets are predicted to have large pt due to the hard 
scattering of the two partons. Observation of this large pr jet production was one 
of early confirmation of QPM in processes other than the DIS. 

Quantitative information about the parton densities can be obtained using jet 
production, however, only in the framework of perturbative QCD. The qualitative 
success of the QPM is formally justified only in the context of perturbative QCD 
(pQCD) any way. Use of jet-production as a probe of parton densities is a little 
complicated as both quarks and gluons in the initial state contribute and over most 
of the range of px neither quarks or gluons dominate. In the framework of pQCD 



one can show that |2^ the large px jet production in hadronic collisions is given by. 



dpTdyidy2 
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f{pT,yi,y2) , 





(4.34) 

where fji = ^in I -g^^pr I is the rapidity of the ith jet in terms of its energy Ei 

and longitudinal momentum PI. The momentum fractions a;i(a;2) of the hadrons 
if 1 (if 2) carried by the partons are related to rapidities yi, 7/2 and pt of the jets and 
f{PT,yi,y2) is a function approximating the dominant part of all the subprocess 
cross-sections da / dp'^{piP2 — ^ P3P4) where Pi stand for partons. Measurement of 
large p^-jet, triple differential cross-section of eq. ( ^4.34 ) by UA-1 collaboration in 



pp collisions |2^ allowed extraction of F^{x) and its comparison with the parton 



densities extracted from DIS, has tested the above picture quite well. 

4.6.2 Drell-Yan process 

Another hard process which was studied in the early days of QPM and does not 
involve pQCD in the subprocess (apart from its role as a justification of QPM) is 
the production of a massive fi'^fj,~ pair with zero px, via qq fusion |Q. The hard 
subprocess, shown in fig. pl|, is 



qq^l* ^ /iV" • (4-35) 

The process at the level of hadrons is 

iiiii2 ^ /x+/i"X . (4.36) 

Let us take the special case when Hi = p, H2 = p ■ Then the total cm energy 
available for the subprocess (say) qP + (f ^ fi~^fi~ is given by, 

^1+,- = S = {P,P + Pgpf = {XlP" + X2P^")2 

~ 2x1X2 P^-P^"~xiX2(P^ + P^)^~xiX25 , (4.37) 

where S is the centre of mass energy of the pp system. The process shown in fig. 
2U| is an electroweak variant of the same. The /^^/i" pair so produced has no net 



momentum transverse to the relative direction of motion of p and p, as the transverse 
momenta of the initial state partons is negligible in the QPM. 

The cross-section for the Drell-Yan (DY) process obviously reflects the quark and 
anti-quark content of both the beam and target. Since the cross-section for fi'^fi~ 
pair production is proportional to e^, the information available from its study is 
essentially the same as that obtained from a study of DIS. In the QPM some simple 
relations exist between DY cross-sections for a fixed target and different beams, e.g., 
the valence quarks in vr"*" are u and d, while those in vr" are u and d. Since the valence 
quarks in proton are u and d, the dominant contribution to DY pair production 
with 71^(71^) beam and proton target comes from d^d'^ [u^u^ ) annihilation. For 



an isoscalar target like Carbon nucleus, we know from our discussion in the earlier 
section, dP = u'~^ . Hence we expect 

a{Ti-C) 4 

This very basic prediction of QPM, which tests the ideas of valence quarks, sea- 
quarks and the quark-content of 7r+/7r^/p etc. as given by the static quark-model, 
was tested quite adequately in the early experiments and thus played an important 
role in establishing the QPM. 

Going outside the framework of plain QPM, in the context of perturbative QCD, 
one can also compute production of massive /x"*"/!" pair, at large pt, in hadron-hadron 
collisions. It takes place via the hard scattering process, 

Qq -^ 7*9 -^ 9l^^f^~] <19 -* Ql* -^ g^+/i". (4.39) 

As discussed for jet production, the final state light parton appears as a jet 
in each case. The 7* (equivalently the ^~^fi~ pair) produced in the hard 2—^2 
scattering subprocess carries a large pt, which is balanced by the jet, even though 
the initial state partons have zero pt- The process at the hadron level is 

Hi + H2^ fi^ + fx^ + jet + X, (4.40) 



where the //"'"/i pair is produced with non-zero pt- It is clear from eq. (|4.39| ) that 



the large pt DY h~^h pair can yield 'direct' information about the gluon content 
of the hadron. 

4.6.3 Direct photon-production 

A related process which can also probe the gluon content of a hadron is the 'direct' 
photon production |2^ via the subprocess 



QQ -^ 91, 19 -^ qi- (4.41) 



The only difference from eq. ( [4. 391 ) is that the photon in the final state is a 'real' 



photon; i.e., the invariant mass of the photon here is zero. The nomenclature 
'direct' (or prompt) has to do with the fact that this photon is produced on the 
short time scale of the hard scattering process as opposed to the photons which 
are produced in the decay of the final state hadrons like 7r° and hence appear on a 
longer (hadronisation) time scale. Of course the latter are produced in much more 
plentiful numbers as the corresponding cross-sections are much higher. Therefore, 
a study of 'direct' or 'prompt' photons is experimentally quite challenging. But the 
effort is well worth it, as one can isolate certain kinematic regions of the final state 
where the cross-section is dominated by the gluon content of the colliding particles 
and hence the measurement offers an almost 'direct' probe of the gluon densities. 



4.6.4 Heavy quark production 

An even better probe of the gluon densities is provided by production of heavy quark 
p^ (charm, bottom etc.) in hadronic collisions via the subprocesses, 

gg^QQ; qq ^ QQ , (4.42) 

where Q stands for c/b quark. This final state has a very distinctive signature and 
can be separated from the very big background of jet production (from the same 
initial states) with comparative ease. The corresponding process at the hadronic 
level is 

H1H2 -^ QQX . (4.43) 

This heavy quark production is dominated by gluon densities due to dynamical 
reasons that make the corresponding cross-sections larger. The situation can be 
further improved (from the point of view the determination of the gluon densities) 
by considering photoproduction of heavy quarks P3 in the process 



-fh-^QQX . (4.44) 

In this case the basic subprocess is 

jg ^ QQX . (4.45) 

This is just an analogue of the corresponding process for the hadronic production 
of heavy flavour, where a gluon is replaced by 7. Thus a study of photoproduction 
(production with incident photons) of heavy quarks can give pretty good information 
about the gluon content of the target. 

A process with smaller cross-section but a more distinctive signature is the pro- 
duction of a bound QQ pair (quarkonium) instead of a free QQ pair. The processes 
which contribute to the hadroproduction (production in hadronic collisions) and 
photoproduction of a quarkonium |]26[| are the same as those in eqs. (|4.42|) and 
( [4.45| ) respectively. 

If we relax the restriction to the 2 — i> 2 subprocess which is inherent to the 
the QPM and allow 2 — >■ 3 subprocesses (which is justified by pQCD) then we can 
describe large pt, photo- and hadro-pro duct ion of quarkonium in the above picture 
{i.e., convolution of parton densities with a subprocess). Some of the subprocesses 
are, 

qq -^ QQg, gg -^ QQg, 75- -^ QQg. (4.46) 

A process which is closely related to the above is the associated production of a 
quarkonium with a photon in hadronic collisions, via the subprocess, e.g., 

gg ^ QQj. (4.47) 

This was recently suggested as a probe of gluon densities. The cross-sections are 
quite a bit smaller than hadroproduction of a quarkonium but, due to the associ- 
ated photon, is much cleaner for detection. It is better than the 'direct' or 'prompt' 



photon in association with jet (eq. ( [4.41|) ), again due to the ease of discrimination 
against background. It also has the advantage of being 'directly' proportional to 
gluon densities for certain spin-parities of quarkonia. In the next section we will dis- 
cuss how these processes can be used to glean information about the gluon densities 
in the nucleus. 

5 EMC effect 

5.1 EMC effect : Experimental situation 

As alluded to in the earlier discussion, because of the high energies involved in 
DIS, it was expected that the nuclear structure function ^2^^ should simply be an 
incoherent sum of the structure functions of the individual nucleons. This picture, 



based on the impulse approximation, depicted in fig. ^ was implicitly assumed in 
all our discussions of the parton model and various tests of the parton model using 
isoscalar nuclei involved using this approximation. The success of these tests, albeit 
qualitative, supports the assumption. However, experiments which set out to test 
this quantitatively met with a surprise ||^. This experiment compared ^iT with 
F2^ . The original experiment studied the EMC ratio {p^^^{x)) defined by 

P^^^(a^) = 7^^ , (5.1) 

A F2^^(x) 

as a function of the scaling variable x for A = Fe. The deviation of this ratio 
from unity is a measure of the failure of the impulse picture. The initial data gave 
pEMC > 1 for a; < 0.3 and a suppression for 0.3 < x < 0.8. Since then a lot of DIS 
experiments verified this non-trivial nuclear dependence of F!^^ {p^^'^'^{x) ^ 1) for 
a wide range of nuclear targets, with a variety of lepton types (e~ beams, /x beams, 
u beams) and over a wide range of q^ (4 < g^ < 200 GeV^)^^. The newer data by 



NMCH probed the effect to very low values of x upto x = 0.0035. Fig. E3| shows a 



collection of some of these data [|T2|, ^ |T9|. The data show the following features : 

1. There is no low x{x < 0.2) enhancement of the nuclear structure functions 
which was present in the original EMC data. 

2. There is a definite 'shadowing' effect, i.e., p^^^^ < 1 for x < 0.05 even at 
large values of q^. The suppression of the nuclear structure function Fg"^ w.r.t. 
the nucleon structure function F2 rises with x for x < 0.2. The 'shadowing' 
means that the DIS cross-section per nucleon is reduced as compared to that 
for a free nucleon, due to the presence of the other nucleons in the nucleus. 
This shadowing effect depends only weakly on q^. 

3. For 0.05—0.1 < x < 0.2 p^^^^{x) goes slightly above unity and then falls below 
1 for 0.3 < X < (0.8 - 0.9). This depletion of ^F^^ w.r.t F^^ is ~ 10 - 15%. 

4. For values of x > 0.8 the EMC ratio goes above unity. 



5. The effect does not show any q^ dependence. 

6. Nor does the effect show strong dependence on the mass number A of the 
target. 

7. An experimental measurement of the ratio R = Fl{x) / (2xFi{x)) does not 
show any appreciable nuclear dependence. 

Due to the intuitive appeal of the impulse approximation, observation of this non- 
trivial nuclear dependence of F2^ was almost unexpected by theorists except for 
a suggestion ||28|| of a possible enhancement of the nuclear sea-densities. For the 
same reason, the observation of the EMC effect also gave rise to a large amount of 
theoretical activity and a large number of models for nuclear structure functions. 
Considering that we do not as yet have a credible model (let alone a theory) even 
of the nucleon structure function, it is clear that all the models proposed to explain 
the nuclear structure functions do so only by giving a recipe to calculate the nuclear 
structure function in terms of that of the nucleon. Different models differ in the 
theoretical ideas about the effect of the nuclear environment on the parton densi- 
ties. All these models of course involve parameters some of which are estimated 
and some are usually fitted to reproduce the observed EMC effect. Since the DIS 
experiments probe only the quark-parton densities directly, it is not surprising that 
all the models agree on the form of the nuclear quark-parton densities. However, 
the nuclear gluon densities which are unconstrained by the DIS data are predictions 
of these different models and generally differ greatly from model to model. The dif- 
ferent models differ radically in the physics phenomenon they invoke to explain the 
EMC effect. Hence, to arrive at the correct theoretical understanding of the EMC 
effect, it is essential to be able to distinguish between the various models. This can 
be done effectively if one can extract nuclear gluon densities and hence measure the 
gluonic EMC ratio 

A gP{x) 

There is yet another reason which makes such a determination imperative. This 
has to do with the signals of the Quark-Gluon-Plasma (QGP) mentioned in the 
earlier sections. To assess the observability of any 'hard' signal of QGP formation 
in Heavy ion collisions it is absolutely essential to understand the contributions to 
the 'hard' final state under question, coming from a combination of the nuclear 
dependence of the parton densities and perturabtive QCD. This of course requires 
a good knowledge of pg defined above. 

The discussions of the last section outlined various hard processes other than the 
DIS which can be used to glean information about the gluonic EMC ratio. Hence, 
a study of correlation between the nontrivial nuclear dependence of the structure 
functions (the EMC effect) and the A-dependence of the different hard processes 
such as large pt jet-production, DY p^p~ pair production(including large pt DY), 
electro- and photo-production of quarkonia as well as their production in hadronic 
collisions etc. can help shed light on the EMC effect. It is worthwhile to ask at this 
point whether there exists any evidence of a nontrivial nuclear dependence for the 



abovementioned hard processes, before we turn to a discussion of such correlations. 
Below we first discuss some of the models that have been suggested to explain the 
EMC effect, summarize the experimental evidence of nontrivial nuclear effects in 
hard processes other than the DIS and then examine the implications of these data 
for various models of the EMC effect. 

5.2 EMC effect: Theoretical Models 

As mentioned in the earlier sub-section, observation of the till then unexpected 
nuclear dependence of the structure functions gave rise to a large variety of modelsQ 
for the EMC effect. Almost all the models address primarily the region 0.2 < x < 
0.8. The small x 'shadowing' region is interesting and has been recently the focus 
of theoretical discussions [^] but will not be discussed here. Broadly speaking the 
models can be divided into different classes: 



5.2.1 Nuclear Physics based models 

Models based on 'conventional' nuclear physics try to explain the depletion of Fr 



A 
2 

in the valence region as being due to the virtual pions present in the nucleus (as 
a result of the nuclear force). The pions can carry a momentum fraction up to 
Mj^/Mp and hence will cause a depletion of quarks in the valence region and also a 
low X enhancement. This idea [^ almost always will give rise to a enhancement of 
the anti-quark content of the nuclear structure function as the pions contain more 
valence anti-quarks. The nuclear structure function is given by 

F2^(x, q') = f dyfM{y)F^{x/y, q') + f dyUy)F^ix/y, q') , (5.3) 

Jx Jx 

where /at and f-,, denote the nucleon and the pion distribution functions in the 
nucleus. Free parameters of this model are the average number of pions in the 
nucleus and the momentum fraction rj carried by the pions. Of course the form 
of the quark-distributions in the free proton and pion, i.e., F^andF^, are also an 
input to the model. The original idea was extended further by a large number of 
authors [Q including also the effect of A's in the nucleus. 

Attempts IP 1|| were made to calculate the two abovementioned parameters in 
nuclear physics framework by using the measured values of the nuclear separation 
energies. There was much excitement initially as the original calculations seemed 
to yield values of t] and the average number of pions required by the best fit to 
the data on EMC effect. However, since then the issue has been revisited by a lot 
of nuclear physicists and there is no clear consensus about the exact size of the 
nuclear binding contribution to the EMC effect; but it is fair to say that it can, at 



the most, account for only 10-15% of the EMC effect ||3^. Since the model always 
predicts an enhancement of the anti-quark content and also in the low x region, the 
experimental information on the nuclear dependence of hard processes other than 
the DIS constrain these types of model rather strongly, as we will see later. 



5.2.2 Cluster models 

Several authors conjectured that the nontrivial nuclear dependence of the structure 
function can be understood in terms of nucleon clusters inside the nucleus. The idea 
is that clusters of nucleons deconfine in the nucleus and share quarks with each other. 
So the quarks belonging to these clusters now occupy the volume of the cluster rather 
than that of a single nucleon. Idea of such clusters was initially suggested to explain 
the deep inelastic data from '^He in the q^ range 1 < g^ < 4GeV^. The probabilities 
of N-quark cluster formation can be computed theoretically. A feature of the cluster 
models is the existence of large momentum partons in the nucleus. For an N-quark 
cluster, the variable x = q^ /2MpV can take values up to A^/3. Thus the fractional 
momentum x carried by the quark can be greater than unity in this case. The cluster 
models (which involve 6 or 9-quark clusters) invoke only a partial deconfinement of 
the quarks in the nucleus. Furthermore the quark distribution functions in the 
free-nucleon and in the N-quark clusters are both input functions. There exist 
also another class of models where the quarks from the nucleons are postulated to 
deconfine to the whole nucleus and these deconfined partons are assumed to form a 
gas. These 'gas' models then compute the parton distributions of the 'deconfined' 
quarks in terms of the parameters of the model. Below two representative models 
of this kind are discussed. 

1) The Gas Model 

The parton densities of a nucleus with atomic number A are defined in the Gas 
model |^3[| as a sum of two components. 

f,/A{x) = (1 - u)h,M{x) + \y. ^'■(1 - ^)^"7r(^; -"' T) . (5.4) 

The first term, occurring with a weight (1 — tu), is that for a free nucleon parton 
density after corrections for its Fermi motion inside the nucleus. The second com- 
ponent is written in terms of thermal distributions of momenta at a temperature T, 
leading to the following functions ff^^X] fi, T): 

(5.5) 
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(5.6) 



Here tq = 1.2 fm, = Mx/2T, z = exp(— /x/T), where M is the nucleon mass and 
Li2{x) is the Euler dilogarithm function. Using the constraint on the total baryon 
number of the nucleus to eliminate the chemical potential /i, one has two model 
parameters, T and u, for each nucleus. 
Using the CDHS parametrisations[p3]. 



(x) = l.l(l + 3.7x)(l-x) 



xap{x) = 0.17(1 -x)^-^^ , (5.7) 

xfg/p{x) = 2.62(1 + 3. 5x)(l-a;)^-^ , 

for -F^x) (here ap{x) is the total sea density) and the data on p{x) = F.^{x)/AF2{x) 
these parameters have been fixed [0^] for many nuclei, including the ones used for 
the E772 experiment. The corresponding Pg{x) is then predicted uniquely. They 
are summarized in Table §. 

It may be mentioned here that, using a different set of structure functions for pro- 
ton instead of eq. (|5.7|) necessitates a time-consuming re-analysis of the EMC-data 
to obtain T and u. For this reason we have not used more recent parametrisations 
for the proton structure function and we are also constrained to use different proton 
structure functions in different models. 

2) The six-quark cluster model 



The six-quark cluster model[^ is also a representative of the two-component models 



for EMC effect. In this model, it is assumed that when two nucleons get closer to 
each other than a certain critical radius they merge together to form a six-quark 
cluster. By assuming the probability to form higher clusters to be negligible, the 
remaining model inputs are the probability of forming such a cluster and the form 
of the parton distributions in the 3-quark and 6-quark clusters. The latter are 
chosen using the quark-counting rules and the constraints of i) normalisation of 
valence densities (an N-quark cluster has N valence quarks), and ii) conservation of 
momentum. It is further assumed that the average momentum carried by the sea 
partons is the same for the three- and six-quark clusters and that it is ~ 0.2 of that 
of the gluons. The forms of nuclear densities per nucleon in this model are, 

fMx) = (1 - e)/,3(x) + '- /,6 (0 , (5.8) 



where e is the probability to find a six quark cluster which increases with A 
and the subscripts denote the cluster size. The values of e used are given in Table 
§. Specific choices]^ for the valence density V{x) = fuvi^) + fdvix)^ sea density 
S{x) and the gluon G{x) for an A^ quark cluster (A^ = 3,6) from Ref.|3^, ^ are 
given by, 

N -1 
xSn{x) = ^p^3^ («iv + 1) (1 - xr- (5.9) 

xfg,N{x) = xGn{x) = ^l^^jf_^l^ (CN + 1) (1 - XY^ , 

with as = 9, a^ = 11, cs = 7, and ce = 10|3^. Here B is the usual Euler Beta 



function and S]\i{x) represents the sum of sea quark densities over all fiavours. With 
a further assumption of fs,N{x) = | fu,N{x) = ^ fd,N{x), the u distribution for 



Table 4: Model parameters for nuclei used in E772 experiment for gas model(T,a;) 
and six quark cluster model(e). 



A 


T (MeV) 


UJ 


e 


12 


54 


0.069 


0.112 


40 


47 


0.057 


0.170 


56 


45 


0.117 


0.186 


184 


42 


0.132 


0.230 



the A^-quark cluster is given by fu,N{x) = | S^^x). The six quark-cluster model 
described is an updated version (in their choice of the parameters e and the input 



proton densities) of the original six quark-bag model [38 



5.2.3 Rescaling models 

A large class of models | 



40 



41| try to model the effect of the nuclear medium 
in terms of the different length scales associated with the nucleon and the nucleus. 
The precise fashion in which it is done varies from model to model. In the g^ 
rescaling models the nuclear parton densities at a scale g^ are obtained from the 
parton densities in a proton at the same g^ by evolving them to a scale ^aq"^, ^-c, 
the nuclear parton density per nucleon fi/A^r, q^) is given by 

fi/A{x,q^) = fi/p{x,U(f) . (5.10) 

Here the results for the rescaled nuclear densities as obtained in Refs. PO] , ^ are 
shown where ^a = ^^"^ where the starting nucleon parton densities are taken [4T| to 
be a parametrisation of the EMC Deuterium data at q^ = 20 GeV^. It should be 
added here that the rescaling models along with the nuclear physics based models 
always tend to enhance the nuclear structure function at small values of x. 

There exist also hybrid models |^ which combine the ideas of both, the rescaling 
and the cluster models. In these types of models the nuclear parton densities(per 
nucleon) are given by, 



fi/A{x,q^) = fi/pix/aA,^Aq'^ 



(5.11) 



The two parameters are introduced to model the change in the scale in the nuclear 
case as well as the the possibility of cluster formation. The two parameters are then 



fitted to reproduce the data on Pemc- Of course tlie fitted values of tlie parameters 
depend on tlie choice of the paranietrisation for the parton densities in the proton. 
The values we obtained [^] are 

a^ = 0.012, a = ^°-'- 

5.3 Comparisons of different model predictions for the gluon 
density 

In fig. ^ are shown the fits to the data obtained in, e.g., the rescaling and the 
hybrid models. This makes the point that the fits to the data on Pemc in different 
models are all of the same quality and all have similar quark parton distributions. 



Fig. |5.4.2| shows the expectations for the gluonic EMC ratio Pg of eq. (|5.2| ), for some 
of the models of the EMC effect discussed above. It should be noted here that the 
different fits to the data use different parametrisations for the proton densities and 
hence it is more meaningful to compare the predictions for pg for the different models 
of the EMC effect rather than the absolute gluon densities. The figure shows clearly 
that the differences in the predictions of the various models are indeed sizable. 

5.4 A dependence of the hard processes 

5.4.1 Experimental situation 

Even before the EMC effect was discovered 0, there existed a few experiments 
which reported an anomalous nuclear enhancement of cross-sections for large-p^ 
particle/jet production |Q and DY p^ p~ pair production |^ with nuclear targets. 



The experiments parametrised the cross-section for nuclear targets with a beam B 
as, 

A ^ ' 

Similarly the ratio of the differential cross-sections, e.g., is parameterised in 

dpT 
terms of a{pT)- Again for no nontrivial nuclear dependence we must have a = 1. A 

deviation of a from unity signals an anomalous nuclear effect. It should be noted 

here that due to the heavy nuclear targets that are used, a small deviation from unity 

for a means a rather large difference between the (per nucleon) cross-section with 

the nuclear and the free nucleon target. The initial experiments [Q reported indeed 

very large values of a increasing with the pt values reaching 1.8 at the highest pt = 6 

GeV. The rise was seen for both p and 7r+ beams. However the jet-like character 

of these data were questionable. In case of the non-jet like data, a large nuclear 

enhancement can also be caused by final state multiscattering effects, which can 

conceivably be larger for nuclear targets as opposed to the nucleon target. Recently 

there has been more data on the nuclear dependence of large pt jet production 

The data show the following features: 



1. The data do show an anomalous nuclear enhancement, i.e., a values bigger 
than unity. The enhancement increases with the transverse energy Et, which 
is essentially a measure of the tranverse momentum pt- 

2. The 'jettier' events give smaller values of a than the non-jet-like events for the 
same value of Et- 

3. The 'jets' seen carry a large fraction of the available cm. energy, y/S. 

The second feature above is consistent with the much larger values of a quoted 
by the earlier experiments. However, in the case of the newer data, due to better 
characterisation of the jet-like nature, multiscattering cannot be invoked to explain 
the anomalous nuclear enhancement and hence has to be interpreted as a reflection 
of the enhancement of Fg^"^ over F2^ in certain x regions. Recently, some data on 
large pt particle production has become available which shows only a modest rise 
of the cross-sections with the mass number A. 

High statistics data on the A dependence of the J/ip and T production |^ 



9| as well as on the large pr DY {fi'^fi~ pair) production [Q have become 
recently available. Experimental data indicate a nuclear suppression of the total 
cross-section. Differential cross-sections are studied in two variables: pt and xp- 
Xp is given by 2Pi/yS, where Pl is the longitudinal momentum of the quarkonium 
or the n^n~ pair and y/S is the total cm. energy, ai^px) and a{xF) show a modest 
Pt and xp dependence respectively. For the large pt DY pairs the a values are 
mostly in the vicinity of unity as opposed to the very large nuclear enhancement 
reported by earlier experiments [^. 

As discussed in the earlier section, yet another interesting probe of the gluon den- 
sities is the photo- and electro-production of quarkonia. The early FNAL and SLAG 
datapTI reported suppression of the lepto-/photo-production cross-sections of the 



J/^p as opposed to the EMC datafS^ which reported an enhancement. The situation 
was clarified by the latest NMC experiments ||53|| where a modest A-dependence of 
the differential distributions of the J/ip production has been reported. 

Thus the experimental studies of the nuclear dependence of the different hard 
processes show quite different behaviours, viz., 

1. Large pt jet production shows considerable nuclear enhancement. 

2. Hadronic quarkonium production shows a suppression for nuclear targets. 

3. The DY n~^n~ pair production and large pt particle production show a very 
modest (almost nil) nuclear dependence. 

From our discussion in the earlier section we know that the dominant subprocesses 
in each of the above hard processes involve different initial state partons. Also 
as a result of the different kinematical conditions the different experiments probe 
parton densities at different values of x. Hence a demand that a given model of 
nuclear parton densities explain all these data on the nuclear dependence of the hard 
processes consistently can indeed constrain these models considerably and help us 
discriminate among them. 



Table 5: Expected values of a in different models compared with the data [@ 





Data 


Gas 


Six-Quark 


Hybrid 


Rescaling 


Config. 


Ref. H 


Ref. g 


Ref. H 


Ref. m 


Ref. [H 


Et > 15 GeV 












P>0.8 


1.14 ±0.02 


1.16 


1.10 


1.03 


0.98 


(Jet-like) 












(with P cut) 












Et > 18 GeV 












Jet-like 


1.45 ±0.01 


1.47 


1.16 


1.06 


0.98 


without P cut 













5.4.2 Comparison of the model predictions with data 

In this section we present some of the model predictions for the A-dependence of 
the different hard scattering processes and compare it with the abovementioned 
data. It should be added here that for consistency, one has to use different proton 
densities while calculating predictions of the different models of the EMC effect for 
the nuclear dependence of various hard scattering processes. 

1) A dependence of jet production 



The first process we consider is large pr particle and jet production with nuclear 
targets. As mentioned in the description of the data, the interesting features of 
the data are that the high-p^ jets seen carry a rather large fraction of the total 
cm. energy, v^- The large value of a implies therefore that the nuclear gluon 
distributions are somehow harder compared with that in a nucleon. This points 
towards cumulative or cluster effects which predict an extension of the nuclear parton 
densities beyond Xbj > 1.0 as opposed to the rescaling models. It can be proved 
on quite general grounds |Q, using simply the experimental information on the 
signature of dp^^^'-^ /dx, that da/dpT will always be less than unless the nuclear 
structure functions extend into the region beyond Xbj = 1, the so called cumulative 
region. These qualitative expectations are indeed borne out by a comparison of the 
model predictions with the data. Table || shows the data along with the predictions 
of different models for the EMC effect. 

A cut on the planarity P helps to choose jet-like events for events with lower 
Ej-- For larger values of E^ the jet-like nature of the events is clear and no such 
cut is required. As one can see, the rescaling model fails to give an enhancement 
of the jet cross-sections. Furthermore, the predicted value of a^ fails even to show 
the A" behaviour. The cluster models [^ |3^, |^ all give values of a bigger than 
1. Stronger the cumulative effects higher are the values of a predicted. Here it is 
worth pointing out that we had used the older version of the six-quark bag model. 
We see that these data already seem to prefer cluster models over the rescaling type 



models of the EMC effect. 

2) Nuclear dependence of the large px quarkonia and /x^/i~ pair pro- 
duction 

Next we discuss a comparison of the data on hadroproduction of large pt quarkonia 
{J/'ip and T) and /i"*"/!" pairs with model predictions. We choose here FNAL data 
||48| , 1^, ^ for its high statistics, although comparisons with the earlier data with 
pion beams ||T7[] do exist [^, ^ . 

The E772 experiment has provided data for the ratio 

dpT I dpT 

with an XiT'-cut of 0.15 < xp < 0.65 on the J/ip^s, while for the T-production cross 
sections, they chose to present only a^pr), where 

dajpA -. TX) ^ ^^^,^^, da{pp-.TX) 
dpT dpT 

with a corresponding xp-cvX for T of —0.2 < xp < 0.6. The nuclei used were 
carbon, calcium, iron and tungsten. One can compute each of the individual pt- 
distributions in eqs. (|5.13| - |5A^ ) incorporating these Xir-cuts. However one needs to 



use a specific model for the quarkonium formation. Fig. |26| exhibits the results for a 
specific model of hadronisation of the quarkonium, for all the four nuclei along with 
the corresponding data from the E772 collaboration. One sees that for the lighter 
nuclei both the two-component models, namely, the gas model and the six-quark 
cluster model, describe the data rather well. For the tungsten nucleus, however, 
none of the models seems to be in agreement with data. 

Fig. |2^ shows a comparison of the model predictions for a values for T-production 
with the E772-data. One sees a similar general agreement for the gas model and 
the six-quark cluster model as for J/'?/'-production at moderate values of pt- At 
the largest pt-, however, the E772-data rise too sharply compared to any model and 
could possibly indicate that these models tuned to earlier large a;-data have to be 
better tuned to perform well in the small x-region. 

The discrepancy at large px values for the T-data and for the tungsten target 
for the J/\E' production do expose the inadequacy of all the three models of the 
EMC effect and the corresponding parametrisation of the nuclear parton densities 
considered here but the general agreement in other cases tells us that the data can 
indeed be described in terms of the structure function effects in general and the data 
are accurate enough to allow discrimination between different models of the nuclear 
structure functions. 

The proton-induced dimuon pair production was studied over a wide range of 
xp and Pt values. The data on the ratio of the integrated dimuon yield for different 
nuclei were compared with theoretical predictions, obtained by using the qq annihi- 
lation process, for various models of the EMC effect. It seemed |BD| to rule out the 



6-quark cluster model[^]. However, a later comparison!^ with an improved ver- 
sion of the model, showed that this model too can be consistent with the information 
on the ratio of the integrated dimuon yields. 

Experimental information [^ is also available for the ratio 

i?^^ = ^^ (M - /^V^x) I ^^{vv^ y^^yrx) , (5.15) 

where da^^ /dpT is the differential DY cross section integrated over the continuum 
region (avoiding the resonances) 4 < M^+^- < 9 GeV and M^+^- > 11 GeV, with 
xf > 0. 



Fig. ^ exhibits the results of a computation for the four different nuclei and the 
three different models with the corresponding data. Again we see that, similar to 
the case of resonance production, the general trends of the data are well described 
by the model predictions for the gas model and the 6-quark cluster model. 

Thus in conclusion we see that already the available experimental information on 
quarkonia, dimuon and large pt jet production indicate that two component models 
of the EMC effect seem to be preferred by the data. Further experiments with direct 



photon production with nuclear targets |5^ or associated production of J/\E' and 



7 |^8|| will help in this direction even more. But what is important to note that it 
is possible now to tune the nuclear gluon densities in different EMC models using 
the data already available. This should go a long way in helping us understand the 
physical origin of the EMC effect as well as help us estimate the QCD backgrounds 
to the 'hard' signals of QGP formation even better. 
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Figure Captions 



Figure 1: Electromagnetic scattering process eA -^ eA along with the four momenta 
assignment for the various particles that are involved. 

Figure 2: Kinematics of the elastic scattering process eA — > eA in the laboratory 
frame. 

Figure 3: Scattering of an electron e~ from the nuclear charge distribution. 

Figure 4: Feynman diagram for the elastic scattering process ep —>■ ep for a pointlike 
proton. 

Figure 5: The form-factor G\.j(q'^)/fip as measured in SLAC experiments as 
a function of g^ (data taken from Ref. Ull)- The solid line is the curve 

(l + gV0.71(GeV2))"'. 

Figure 6: Kinematics of the inelastic scattering process ep — > eX. 

Figure 7: Allowed region in the q'^ — u plane for the elastic, quasi-elastic and inclusive 
inelastic scattering. 



Figure 8: W2 = M^ 1^2(1^, Q'^)/iTa as a function of g^and u (Figure taken from Ref. 



|T6[|). The details of the source of the data etc. to be found in Ref. |jT6 



Figure 9: Data on VW2' as a function of q^ for u = 1/x = 1/4 (taken from Ref. |0 
Different data correspond to different scattering angles as indicated in the figure. 



Figure 10: Parton model picture of the DIS scattering of an e off a proton target. 



Figure 11: DIS scattering of a i/ off a proton target. 



Figure 12: Schematic drawing of F2 as a function of x for different values of q^ in 
the range 0.01 - 200 GeV^ (taken from Ref. g 



Figure 13: Data on e — a scattering at g^ = 0.08 GeV^ (taken from Ref. [Q)- A is 
the elastic peak for a particle while the elastic proton peak is shown by the dashed 
line which corresponds to x = q^ /2ma v — 0.25. The portion BCDE indicates the 
momentum distribution of the nucleons in the a particle. 

Figure 14: Parton model. 

Figure 15: Impossibility of the backward scattering for the helicity preserving inter- 
action. 

Figure 16: Data on -pwrfy as a function of x (data taken from Ref. p2|). 

Figure 17: The ratio R of eq. ( |4.23| ) as a function of g^ (taken from Ref. 



Figure 18: Ratio of F2^^ from the CCFR and CDHSW experiment to F2 from 
different fi DIS experiments. The fi structure functions are normalised so thata 
ratio of unity means mean square charge of -j^ (figure taken from Ref. 1^). See 



Ref. pOf for more details. 

Figure 19: Data on F2^ — ¥2^ from NMC as a function of x (data from Ref. [p^ ) 

Figure 20: Parton model picture of the W^ production in pp collisions. 



Figure 21: DY process of production of a fi~^fi pair production in hadronic colli- 
sions. 



Figure 22: Impulse approximation for the nuclear DIS scattering. 



Figure 23: Compilation of the data on EMC effect (from Ref. [|T2|)- Details of the 
data available in Ref. [^ 



Figure 24: Fits to the data on pemc in some models of the EMC effect [^ . 



Figure 25: The ratio pg of eq. ^^ as a function of x for some models of the EMC 
effect m. 



Figure 26: E772 data (inverted filled triangles) on the ratio R'^/^ of Eq. p.13 com- 
pared with the predictions for the gas model (squares), six-quark cluster model 
(circles) and the rescaling model (open triangles) of the EMC effect [^ . 



Figure 27: E772 data on a{pT) of eq. |5.14] compared with predictions of the three 
different models of the EMC effect mentioned in the text. Notation is same as in 
Fig. 
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Figure 28: E772 data on the ratio R^^ of Eq. |5.15| compared with predictions of 
the three models of the EMC effect. Notation is same as in Fig. El 



A Notations used 

If A and B are two four vectors given hj A = (a, itto), B = {b, iho) then, 

A ■ B = a ■ b — aobo = a ■ b + 0464. 
The 7 matrix algebra is given as 

7^ = 7/. 
7^ = 1 
1 

75 = 71727374 = -^e^^,^pa1^,1ulpla 

{75,7.} = ; /i= 1,2,3,4 

75 =75 

(^/.u = ^ [7m' 7i^] = -^ 7^7;^ (/^ 7^ J^) 

The Dirac equation for a particle of mass m is given in this metric by 

(7^9^ + m)^ = 0. 

In momentum space this becomes: 

{i f + m)u{P) = where ^ = j^P^ 

P stands for the four-momentum of the particle and u{P) is the free particle spinor. 
In this metric we have 

Y^ Ua{P, s)Ui3{P, S) = {-if + m)ap. 
s 

The trace theorems in this metric are given by: 

Tr(7^7^) = 4 6^y 
Tr [7M7!^7p7a] = 4 [5f,Jp„ + 5^„5yp - 5pp5^,„] 
Tr [7M7i^7p7a75] = 4 e^^p^ 
Tr [71 ■ ■ ■ 7„] = (n odd) 
Tr [757^ ■ ■ ■] = if 7^ • ■ ■ has less than 4 7 matrices 

One more relation required in calculation of (|7Wp) is 

74[(7 ■ ^)(7 • 5). ...(7 • L)]\, = (-7.L) • ■ ■ (-7 ■ A); 
where A,B---L stand for four vectors representing four momenta. 



B Problems 

1 What is the energy of the probing e~ beam that will be required to probe the 
structure at a distance scale (say) 10^^'^ m? 

2 Calculate the form factors F[Q'^) for the following charge distributions: 1) 
p{R) = S%R), 2) p{R) = ^exp{-mR)/R and 3)^exp{-mR) where m > 0. 

3 Show by explicit calculation that the factor cos^(6'/2) in eq. ( |1.42| ) corresponds 
to the helicity non-flip amplitude. 

4 Calculate the expression for the differential cross-section ^ for the electro- 
magnetic scattering of an electron incident on a target which is 

1) spin 0, pointlike particle 

2) spin 1/2, pointlike particle 

3) spin 1/2 particle which is not pointlike. 



The expression for the cross-section is given by eq. (|2.4| ) The matrix element 



is given by eq. (|2.3D where the electromagnetic current of the electron is given 



by eq. ( |2.9| ), and the electromagnetic current of the target in cases 2 and 3 
are given by eqs. ([2.9[) and (|2.14D respectively, whereas the electromagnetic 



current for the case 1) of a spinless, pointlike proton is given by 

J', = ^qpiP4 + Pi), 

where all the notations are as given in the lecture. The answers are given by 
eqs. (U),(|3D and (|2T7| ). 



5 Show that considerations of gauge invariance and Lorentz invariance restrict 
the form of H" as given by eq. ( |3.45|) , by following steps analogous to those 



used in deriving eq. ( ^.221 ). Using this, show that the expression for the 
differential cross-section is given by eq. ( p.46| ). 

6 Calculate the expressions for F2^ and F^^ given that, 

da , GpS 

ay IT 



Assume that the Callan-Gross relation given by eq. ( |4.5| ) is satisfied for the 
Neutrino structure functions F2^ and Fi^ and use the expression for the 
double differential cross-section for the inelastic up scattering given by eq. 



7 Show that in parton model, for an isoscalar nuclear target 



FT 



1 



4:K{x) 



Uy{x) + dy{x) + 4 K{x) 



where K{x) is the SU(3) symmetric sea densities defined in eq. (|4.15|) . 
8 Derive eq. ( P5| ). 



